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Abstract. This review article discusses limit distributions and vari- 
ance bounds for particle current in several dynamical stochastic systems of 
particles on the one-dimensional integer lattice: independent particles, in- 
dependent particles in a random environment, the random average process, 
the asymmetric simple exclusion process, and a class of totally asymmetric 
zero range processes. The first three models possess linear macroscopic flux 
functions and lie in the Edwards- Wilkinson universality class with scaling 
exponent 1/4 for current fluctuations. For these we prove Gaussian limits 
for the current process. The latter two systems belong to the Kardar-Parisi- 
Zhang class. For these we prove the scaling exponent 1/3 in the form of 
upper and lower variance bounds. 
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1 Introduction 

This review article investigates the process of particle current in several con- 
servative stochastic systems of particles that live on the one-dimensional 
integer lattice. In conservative systems particles are neither created nor de- 
stroyed. On the macrosopic, deterministic scale the particle density p(t, x) 
of such systems is governed by a partial differential equation of scalar con- 
servation law type: pt + H(p) x = 0. The results covered in this article 
concern the fluctuation behavior of the current, and include some precise 
limit results and some coarser order-of-magnitudc bounds. Some back- 
ground on advanced probability theory is assumed. For readers with little 
probability background we can also suggest review article [Scp08] where an 
attempt was made to explain some of this same material for a reader with 
background in analysis rather than probability. 

The particle processes studied fall into two categories. We can define 
these two categories by the slope of the Eux function H: (i) processes with 
a linear flux and (ii) those with a strictly concave flux. By definition, the 
flux H (p) is the mean rate at which particle mass moves past a fixed point in 
space when the system is stationary in both space and time and has overall 
density p of particles. The processes we study are asymmetric or driven 
in the sense that the particles have a drift, that is, a preferred average 
direction. This assumption is not necessary for the results for systems with 
linear flux, but for those with nonlinear flux it is crucial. 

In the statistical physics terminology of surface growth, we can also la- 
bel these two classes as (i) the EW (Edwards- Wilkinson) and (ii) the KPZ 
(Kardar-Parisi-Zhang) universality classes [BS]. Surface growth may seem 
at first a separate topic from particle systems. But in one dimension conser- 
vative particle systems can be equivalently formulated as interface models. 
The connection goes by way of regarding the particle occupation numbers 
rji as increments, or discrete gradients, of the interface height function: 
hi = rji — r]i-i. Then any movement of particles in the conservative particle 
system can be equivalently described as deposition or removal of particles 
from the growing surface. In particular, the current process then maps 
directly into the height function. 

In the EW class limiting current/height fluctuations are described by the 
linear stochastic heat equation Z t = vZ xx + W where subscripts are partial 
derivatives and W represents space-time white noise. In the microscopic 
model the order of magnitude of current fluctuations is n 1 / 4 in terms of a 
scaling parameter n (n /• oo) that gives both the space and time scale. The 
results we give describe Gaussian limit distributions of the current process. 
The examples of systems with linear flux we cover are independent random 
walks, independent random walks in a random environment (RWRE), and 
the random average process. 

In the RWRE case we in particular wish to understand how the envi- 
ronment affects the fluctuations. We find a two-level fluctuation picture. 
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On the diffusive n 1 / 2 scale the quenched mean of the current converges to 
a Brownian motion. Around the quenched mean, on the n 1 / 4 scale, the 
current converges to the Gaussian processes that arose in the case of inde- 
pendent classical random walks, but with an additional Brownian random 
spatial shift. 

The systems with concave flux that we discuss are the asymmetric simple 
exclusion process (ASEP) and a class of zero range processes. In this case 
the assumption of asymmetry is necessary. We consider only stationary 
systems (or small perturbations thereof) , and instead of distributional lim- 
its we give only bounds on the variance of the current and on the moments 
of a second class particle. The current fluctuations are now of order i 1 / 3 . 
This order of magnitude goes together with superdiffusivity of the second 
class particle whose fluctuations are of order t 2 ^ 3 . Here t is the time pa- 
rameter of the process. A separate scaling parameter is not needed since 
we have no process level result. These systems are in the KPZ class. 

We cannot cover complete proofs for all results to keep this article at a 
reasonable length. The most important result, namely the moment bounds 
for the second class particle in ASEP, is proved in full detail, assuming 
some basic facts about ASEP. 

In some sense interface height in the KPZ class should be described by 
the KPZ equation 

ht = vh xx - \{h x ) 2 + W. (1.1) 

However, giving mathematical meaning to this equation has been done 
only indirectly. A formal Hopf-Cole transformation Z = exp(— A^ _1 /i) 
converts (1.1) into a stochastic heat equation with multiplicative noise: 
Z t = vZ xx + \v~ 1 ZW . The solution of this latter equation is well-defined, 
and can furthermore be obtained as a limit of an appropriately scaled height 
function of a weakly asymmetric simple exclusion process [BG97]. Via this 
link the scaling exponents have recently been verified for the Hopf-Cole 
solution of the KPZ equation [BQS]. 

The most glaring omission of this article is that we do not treat the Tracy- 
Widom fluctuations of KPZ systems such as TASEP (totally asymmetric 
simple exclusion process), ASEP or the PNG (polynuclear growth) model. 
In the ASEP section we only state briefly the Ferrari-Spohn theorem [FS06] 
on the limit distribution of the current across a characteristic in stationary 
TASEP. This area is advancing rapidly and requires a review article of 
its own if any degree of detail is to be covered. Tracy and Widom have 
recently written a series of papers where the fluctuation limits originally 
proved for TASEP by Johansson [JohOO] have been extended to ASEP 
[TW08a, TW08b, TW09a, TW09b, TW10]. Another line of recent work 
extends the TASEP limits to more general space-time points and initial 
distributions [BFP, CFP]. 



A notational convention. Z + = {0, 1,2,3,...} and N = {1, 2,3,...}. 
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2 Independent particles executing classical 
random walks 

2.1 Model and results 

Fix a probability distribution p(x) on Z. Let us assume for simplicity that 
p(x) has finite range, that is, 

the set {ieZ: p(x) > 0} is finite. (2.1) 

For each site x £ Z and index k £ N let Xf ' be a discrete-time random 
walk with initial point Xq ' k = x and transition probability 

P(X:*=y\X*' k = z)=p(y-z) 

for times s £ Z + and space points z, y £ Z. The walks {X?' : x £ Z, k > 1} 
arc independent of each other. Let 

v = xp(x) and a\ = ' S ^ j {x — v) 2 p(x) 

be the mean speed and variance of the walks. 

At time we start a random number n x (0) of particles at site x. The 
assumption is that the initial occupation variables 77(0) = {^(0) : x £ Z} 
are i.i.d. (independent and identically distributed) with finite mean and 
variance 

Ho = E[n x {Q)] and a% = Varfo^O)]. (2.2) 

Furthermore, the variables {^(0)} and the walks {X*' k } arc independent 
of each other. If the locations of individual particles arc not of interest 
but only the overall particle distribution, the particle configuration at time 
s £ N is described by the occupation variables ry(s) = {r] x (s) : x £ Z} 
defined as 

y£Z fc=l 



6 



Timo Seppalainen 



A basic fact is that if the initial occupation variables are i.i.d. Poisson 
with common mean (i, then so are the occupation {rj x (s) : x £ Z} at any 
fixed time s £ Z + . This invariance follows easily from basic properties of 
Poisson distributions: the counts of particles that move from x to y are 
independent Poisson variables across all pairs (x,y). 

Let E^- denote expectation under the stationary process with mean \x 
occupations. The flux function is the expected rate of flow in the stationary 
process: 



H(n) = -E M [net number of particles that jump across edge (0, 1) 
left to right in one time step] 

no(x) vo(. x ) 



VfX. 



E E !{*r fc > i} - E E W < °} 

x<0 k=l x>0 k=l 



(2.3) 



The linear flux puts this system in the class where we expect n 1 / 4 mag- 
nitude current fluctuations. Next we define the current process and the 
Gaussian process that describes the limit. 

Let n £ N denote a scaling parameter that eventually goes to oo. For 
(macroscopic) times t £ R+ and a spatial variable r £ R, let 



with 



and 



Y n (t,r) =Y nil (t,r) -Y n , 2 (t,r) (2.4) 
Y nA (t, r) = EE 1 { Jr L»tj > ^ + r ^ > ( 2 - 5 ) 

x<0 k=l 

Y nt2 (t,r) = E E < W (2.6) 

a:>0 fc=l 

Variable F„(t, r) represents the net left-to- right current of particles seen by 
an observer who starts at the origin and reaches point [^^J + L r v^J m 
time \ nt\ . Its mean is 

EY n {t,r) = n E(X lntl - [ntv\ - [ry/ri\) = -ftory/n + 0(l). 

Define the centered and appropriately scaled process by 

Y n (t,r) =n- 1 / 4 (Y n (t,r) -EY n (t,r)). 

The goal is to prove a limit for the joint distributions of these random 
variables. We will not tackle process-level convergence. But let us point out 
that there is a natural path space D 2 of functions of two parameters (t, r) 
that contains the paths of the processes Y n . Elements of D 2 are continuous 
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from above with limits from below in a suitable way, and there is a metric 
that generalizes the standard Skorohod topology of the usual D-space of 
cadlag paths. (See [BW71, Kum08].) 
Let 

1 { x 2 \ f x 

= ^pHv 2 ° XP V 2V 1 ) and ® v2 t x ) = J <Pv 2 (y) d V ( 2 - 7 ) 

denote the centered Gaussian density with variance v 2 and its distribution 
function. Let W be a two-parameter Brownian motion on R_|_ x R and B a 
two-sided one-parameter Browian motion on R. W and B are independent. 
Define the process Z by 

Z(t,r)=^/jMi (p a 2 (t ^ s) (r - x)dW(s,x) 

JJ[0,t]xR ^ g) 



+ °o / <p a 2 t (r - x)B(x)dx. 
Jr 

{Z(t,r) : t £ R+,r € R} is a mean zero Gaussian process. Its covariancc 
can be expressed as follows: with 

^> V 2 (x) = v 2 Vv 2 (x) - x(l - $,2 (x)) (2.9) 

define 

rj ((a, g), (t, r)) - * CTf(t+s) (r - q) - * CT 2 |t _ s| (r - q) (2.10) 

and 

r 2 ((s, g), (t, r)) = * CTfs (-g) + * CT?t (r) - * CTf(t+s) (r - g). (2.11) 

Then 

E[Z(s, g)Z(i, r)] = Mori ((a, g), (i, r)) + a 2 r 2 (( s , g), (t, r)) . (2.12) 

Remark 2.1. Some comments on Gaussian processes defined above. Com- 
parison of (2.8) and (2.12) shows that Ti is the covariance of the dynamical 
fluctuations represented by the space-time white noise dW^-integral, and r 2 
is the covariance of the contribution of the initial fluctuations represented 
by the Brownian motion aoB(-). 
Ti has this alternative formula: 

r 1 ((a,g) 1 (t,r))=i -jL= exp{l(r - q) 2 } dv. (2.13) 



To verify that the process 



C(t,r)= ip a 2 {t _ s) (r - x)dW(s,x) 

' '[0,t]xR 
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has covariance 

EC(fl,g)C(*,r)=ri((«,g),(t,r)) (2.14) 

is straightforward from the general property that for /, g £ L 2 (R d ) the 
white-noise integrals / / dW and J g dW on R d are by definition mean 
zero Gaussian random variables that satisfy 



E 



fdW)( / gdW 



f(x)g(x) dx. 



To say that B(x) is a (standard) two-sided Brownian motion means that 
we take two independent standard Brownian motions B\ and Bi and set 



Bx(x), x>0 
B 2 (-x), x<0. 



B(x) 

To show that the process 

r ) = Vvftir ~ x)B(x) dx 

is a mean-zero Gaussian process with covariance 

ms,qW,r)=T 2 ((s,q),(t,r)) 

this formula is useful: if /, g are absolutely continuous functions on R + 
such that xf{x) — ¥ and xg(x) — > as x —> 00, then 

/>oo 

f'(x)g'(y)(xAy)dxdy= / f{x)g{x)dx. 



It is also the case that the process Z(t, r) is the unique weak solution of 
the following initial value problem for a linear stochastic heat equation on 
R+ x M: 

Z t = ^Z rr + ^W, Z(0,r)=a B(r). 

(Above, subscript means partial derivative.) A weak solution of this equa- 
tion is defined by the requirement 

<{f)Z(t, r) dr — Co / 4>(r)B(r)dr 



4>"{r)Z{s,r)drds + Vm / / 4>(r)dW(s,r) 

[0,t]xK JJ[0,i]xl 

for all 4> £ C^°(M) (compactly supported, infinitely diffcrcntiablc). See the 
lecture notes of Walsh [Wal86] . 

We can now state the main result. 
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Theorem 2.2. Assume the initial occupation variables are i.i.d. with finite 
mean and variance as in (2.2). Then as n — > oo, the finite- dimensional 
distributions of the process {Y n (t,r) : (t,r) £ R + x M.} converge weakly 
to the finite- dimensional distributions of the mean zero Gaussian process 
{Z(t,r) : (t,r) £ R+ x R}. 

The statement means that for any space-time points (ti, r\ ),•••, (ifcjffc), 
this weak convergence of K fc -valued random vectors holds: 

(F n (ti,n), . . . ,Y n (t k , r k )) A (Z(ti,n), . . . , Z{t k , r k )). 

Under additional moment assumptions process level convergence in the 
space D<i can be proved (see [Kum08]). We state a corollary for the spe- 
cial case of the stationary occupation process rj{t). Its proof comes from 
simplifying expression (2.12) for the covariance. 

Corollary 2.1. Suppose the process is stationary so that {r] x (t) : x £ Z} 
are i.i.d. Poisson with mean /io for each fixed t. Then at r = the limit 
process Z has covariance 

EZ(s, 0)Z(t, 0) = ^1(^ + Vt- VW^l) ■ 
v 2n 

In other words, process Z(-,0) is fractional Brownian motion with Hurst 
parameter 1/4. 

2.2 Sketch of proof 

We turn to discuss the proof of Theorem 2.2. Independent walks allow us 
to compute everything in a straightforward manner. Fix some TV £ N, time 
points < t\ < £2 < ■ ■ • < tN S K+, space points r\,V2, ■ ■ ■ ,Tn G K and 
an ^-vector = (6±, . . . , 6jy) £ M. N . Form the linear combinations 

N N 

Yn(0) =Ys 9 i Y n( t i> r i) and Z W = Y,^ 2 ^^' 

i=l i=l 

The goal is now to prove 
Proposition 2.1. 

E[exp{iY n {8)}] ^E[cxp{iZ(0)}]. (2.15) 



Since the random walks and initial occupation variables are independent, 
we can write Y n (6) as a sum of independent random variables and take 
advantage of standard central limit theorems from the literature. 



Y n (8) = n-?J26 i {Y n (t i ,r i )-EY n (t i ,r i )}=W n = ^ u(m) (2.16) 
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with 

N 

u(m)= Y, 9 i( U m(U,ri)l{m<0}-V m (U,ri)l{m>0}), (2.17) 



(2.18) 



and 

U m (t,r) = n~i J2 HXIS? > \ntv\+r^} 

3 = 1 

-n-i{M>P{X™ t > [ntv\ + rVE), 

»7m(0) 

V m (t,r) = n-i J2 ^ X "t 3 ^ [ntvj+rV^j 
j'=i 

-Tr*HoP(X™ t < [ntv\ + ry/n). 

The variables {u(m)} me i are independent because initial occupation vari- 
ables and walks arc independent. 

Let a(n) / oo be a sequence that will be determined precisely in the 
proof. As the first step we observe that the terms \m\ > a(n)^/n can be 
discarded from (2.16). Define 

w* n = < m )- ( 2 - 19 ) 

\m\<a(n)*/n 

The lemma below is proved by calculating moments of the random walks. 

Lemma 2.1. E\W n — W*\ 2 — >• as n -> oo. 

The limit Z(ff) in our goal (2.15) has Af(O,a{0) 2 ) distribution with vari- 
ance 

a(6) 2 = Yl O l d^ T 1 ((U 1 r l ) 1 (t ] ,r 1 ))+a 2 T 2 ((U,r x ),(t J ,r ] )) . 

l<i,j<N 

(2.20) 

The two T-terms, defined earlier in (2.10) and (2.11), have the following 
expressions in terms of a standard 1-dimcnsional Brownian motion B t : 



ri((a,g),(t,r))= / (p[B ajs < q - x]P[B ajt > r - x] 

J — OO 

- P[B a 2 s <q-x, B a i t > r-x])dx 



(2.21) 



and 



T 2 ((s,q),(t,r)) = / P[B als >q-x]P[B aft >r-x]dx 

J —OO 

+ / P[B a 2 <q-x]P[B a 2 t <r-x\dx. 
Jo 



(2.22) 
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Remark 2.3. Turning formulas (2.21)-(2.22) into (2.10)-(2.11) involves 
calculus and these properties of Gaussians: (d/dx)^^ (x) = — $„2(— x), 
$„2(a;) = 1 - $„a(— x) and 



/oo />r 
$ Q 2(x)$^2(r-x)dx = / $ Q 2 + ^(x) 
-oo J —oo 



dx. 



By Lemma 2.1, the desired limit (2.15) follows from showing 

£(e w n)_> e — W 2 /2. (2.23) 

This will be achieved by the Lindeberg-Feller theorem. 

Theorem 2.4 (Lindcberg-Fellcr). For each n, suppose {X n ,j : 1 < j < 
J(n)} are independent, mean-zero, square-integrable random variables and 
let S n — X n ^i + ■ ■ • + X n jt n \. Assume that 

J(n) 

lim J2 E ( x h) = ° 2 



n—*oc 

3=1 



and for each e > 0. 



J(n) 

lim J2 E ( X n jH\X nd \>e})=0. 



n— J-oo 



T/iera as n —+ oo, S n converges in distribution to a Af(Q, a 2 ) -distributed 
Gaussian random variable. In terms of probabilities, the conclusion is that 
for all s e R, 

lim P{S n < s} = - 1 / e- x ' 2/2a2 dx. 

™^°° V27TCT 2 J-oo 

In terms of characteristic functions, the conclusion is that for all t£R, 

lim E(e lts ") = e^ 2 ' 2 / 2 . 

Now to prove (2.23) the task is to verify the conditions of the Lindeberg- 
Feller theorem: 

5^ E{u{m) 2 ) -> <r(0) 2 (2.24) 

|m|<a(n)y / n 

and 

E(\u(m)\ 2 l{\u(m)\ > e}) -> 0. (2.25) 

We begin with the negligibility condition (2.25). This will determine 
a(n). 
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Lemma 2.2. Under assumption (2.2), 

lim V E( |u(m)| 2 l{|u(m)| > s}) = 0. (2.26) 

n—too ' J 

\m\<a(n)*s/n 

Proof. Since 

\u(m)\ <Cn- 1 ^(r hl (0)+^) 
so for a different e\ > and by shift-invariance, 

E E[\u(m)\ 2 l{\u(m)\ > e}] < Ca(n)E[( Vo (0) + Mo) 2 l{?7o(0) > n 1 / 4 

By the moment assumption (2.2) this last expression — > for every e\ > 
if a(n) /* oo slowly enough, for example 

a(n) = (E[(m(0) +mo) 2 1{%(0) > n 1 / 8 }])" 1 '* □ 
We turn to checking (2.24). 

E ^[«(«) a ] 

|m| <a(?i)-\/n 

= E ^ E 1 {m<o}£ ; (C / m(i l ,r-j)[/ m (^,r i )) (2.27) 

l<i,j<N |m|<o(n)Vn 

T 1 {m>0} 

To the expectations we apply this formula for the covariance of two random 
sums: with {Z{\ i.i.d. and independent of a random K G Z+, 

if AT 

Cov(E/(^),E^)) 

»=1 3=1 

= £A Cov(f(Z),g(Z)) + Var(A') S 5 (Z). 
For the first expectation on the right in (2.27): 
E(U m (s,q)U m (t,r)) 

y»)m(0) 

= ,i- 1 /2 C ov E X W > M +<Zv^}, 

^ 3 = 1 

Tlm(0) 

E lffi j >N+rA} 

= n" 1/2 Mo > M + gV^, ^ > [n*«J + r V^) 

- P(XZ > [nsv\ + qy/Z)P(X% > [ntv\ + ryfa) 
+ n- 1 ' 2 olP{XZ > [nsv] + qy/E)P(X^ t > \ntv\+r^i). 



(2.28) 
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Do the same for the V-terms. After some rearranging of the probabilities, 
we arrive at 

£ E[u(m) 2 ] (2.29) 

\m\<a(n)y^n 

^n- 1 ' 2 £ W £ { P ( X nU<l nt iv\+nV^) (2.30) 

\m\<a(n)y/n 

xP(X™ > [ntjv} +r,V^) 



l<i,j<N 



+ ct 2 £ p ( x ^ > L**<«J + ^)PpC. > Ln*i«J + ^) 

— a(n)i/n< m < 
0<m<a(n)^/n 

The terms above have been arranged so that the sums match up with 
the integrals in (2.20)-(2.22). Limit (2.24) now follows because each sum 
converges to the corresponding integral. To illustrate with the last term, 
the convergence needed is 



-1/2 

0<m<a(n)y^n 

n -i/ 2 y rX,.-Ln^ < m | 

0<m<a(n)Vn (2.31) 

n f X nt - [ntv\ m ■} 

xP jk J -<r- — 

\ n \/n J 



/ P[B als <q-x]P[B a 2 t <r-x}dx. 
Jo 



This follows from the CLT, a Riemann sum type argument and some esti- 
mation. We skip the details. With this we consider Theorem 2.2 proved. 
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3 Independent particles in a random envi- 
ronment 

In this chapter we generalize the results of Section 2 to particles in a random 
environment, with the purpose of seeing how the environment influences the 
outcome. In a fixed environment, that is, conditional on the environment, 
the particles evolve independently. But under the joint distribution of the 
walks and the environment, the particles are no longer independent because 
their evolution gives information about the environment. The environment 
is static, which means that it is fixed in time. 

3.1 Model and results 

We formulate the standard one-dimensional nearest-neighbor random walk 
in random environment (RWRE) model and then put many particles in 
a fixed environment. We describe the (known) law of large numbers and 
central limit theorem of the walk itself, and then the (newer) results on 
current fluctuations for many particles. 

The space of environments is f2 = [0, For an environment u = 
{^x}x& G O let {X™' l } m ^ be a family of Markov chains on Z with distri- 
bution determined by the following properties: 

1. {X m ' I }mez,iGN are independent under the measure P u . 

2. P^iX™- 1 = m) = 1, for all m G Z and i G N. 

3. Each walk obeys these transition probabilities: 

PUK+i = x+iixp* = x) = i-PUK+i = W- 1 = *) = 

A system of random walks in a random environment may then be con- 
structed by first choosing an environment u> according to a probability dis- 
tribution P on Q and then constructing the system of random walks {X m ' 8 } 
as described above. The distribution P u of the random walks given the en- 
vironment ui is called the quenched distribution. The averaged distribution 
P (also called annealed) is obtained by averaging the quenched law over 
all environments: P(-) = P ul (-)P(du>). Expectations with respect to the 
measures P, P u and P are denoted by Ep, E u , and E, respectively, and 
variances with respect to the measure P u will be denoted by Var w . We 
make the following assumptions on the environment. 

Assumption 1. The distribution P on environments is i.i.d. and uni- 
formly elliptic. That is, {uj x } x ^z are i.i.d. under the measure P, and 
there exists a k > such that P(lo x g [k, 1 — k]) = 1. Furthermore, 
Ep(Pq +£ ") < 1 f° r some £o > 0, where p x = (1 — u) x )/oj x . 
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These assumptions put the RWRE in the regime where it has transience 
to +00 with a strictly positive speed and also satisfies a CLT with an 
environment-dependent centering. We summarize these results here. De- 
fine a shift map on environments by (6 x uj) y = tu x +y Let T\ = inf{n > : 
X n = 1} be the first hitting time of site 1 <G Z by a RWRE started at the 
origin, and define 

[ntvpj-l 

Z nt (Lu)= Vp (^fl«a,(Ti)-ETi). (3.1) 

i=0 

The asymptotic speed vp is defined in the first statement of the next the- 
orem where we summarize some known basic facts about RWRE. 

Theorem 3.1 ([Sol75, Pet08, Zci04]). Under the assumptions made above 
we have these conclusions. 

1. The RWRE satisfies a law of large numbers with positive speed. That 
is, 

km = — r = vp > 0, F-a.s. (3.2) 

n-s-oo n l + hp{p ) 

2. The RWRE satisfies a quenched functional central limit theorem with 
an environment- dependent centering. Let 

nn/.\ X nt ~ ntvp + Z nt (uj) 9 3r ,,, 

B (t) = -= — , where erf = V P Ep{Va,r ul T 1 ). 

Then, for P-a.e. environment ui, under the quenched measure P^, 
B n (-) converges weakly to standard Brownian motion as n — > 00. 

3. Let 

= where a \ = ^ Var(£ w T x ). 

02 V n 

Then, under the measure P on environments, C™(') converges weakly 
to standard Brownian motion as n — > 00. 

4- The RWRE satisfies an averaged central limit theorem. Let 

m nf.\ X nt -ntv P 2 2,2 

B (t) = -= , where a = a x + a 2 - 

Then, under the averaged measure P, B™(-) converges weakly to stan- 
dard Brownian motion. 

The requirement that Ep(pfy < 1 cannot be relaxed in order for the CLT 
to hold [KKS75, PZ09]. Centering by ntv P - Z nt {uj) in the quenched CLT 
is the same as centering by the quenched mean on account of this bound: 

km p{uj : sup \E u (X k ) - kv P + Z k {uj)\ > = 0, Ve > 0. (3.3) 

n->-oo I k<n J 
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But Z nt (uj) is more convenient because it is a sum of stationary, ergodic 
random variables. 

These properties of the walk are sufficient for describing the current fluc- 
tuations. Assumptions on the initial occupation variables r)(0) = {^(O)} 
are similar to those in the previous section. We will allow the distribu- 
tion of 77(0) to depend on the environment (in a measurable way), and we 
assume a certain stationarity. 

Assumption 2. Given the environment to, variables {r] x (0)} are inde- 
pendent and independent of the random walks. The conditional distribu- 
tion of 77^(0) given u is denoted by P u (^(O) = fc), and these measurable 
functions of w satisfy Puj(n x (0) = k) = Pe^uiivoi^) = k). Also, for some 
£0 > 0, 

E P [EUriA0)) 2+eo + Var^O)) 2 ^ ] < 00. (3.4) 

Let 

Mo = E P [E u (r] x (0))} = E[r? x (0)] and <t 2 = E P [Var w fe(0))] . 

The current is defined as before: 

rim (0) 

Y n(t,r) = V V l{X^ k > [ntv P \ + r^} 



m<0 k=l 



ri m (0) 

Y, E l{X^ k <[ntv P \+rV^}. 

m>0 k=l 



(3.5) 



Now for the results, beginning with the quenched mean of the current. 
This turns out to essentially follow the correction Z nt (uj) of the quenched 
CLT, which is of order ^pri. Assumptions 1 and 2 are in force for all the 
results that follow. 

Theorem 3.2. For any e > 0, < R, T < 00, 

lim Pico : sup \E u Y n (t, r) + fi r^/n + (XoZ-nt M I > eVn f = 0. (3.6) 

n^oo L te[0,T] J 

re[-R,R] 

Consequently the two-parameter process {n- l l 2 E u Y n {t,r) :teR+, ret} 
converges weakly to {— /^o?" + no<j2W(t) : t € R+, r£ 1} where W(-) is a 
standard Brownian motion. 

Next we center the current at its quenched mean by defining 

V n (t,r) =Y n (t,r)-E u [Y n (t,r)]. 

The fluctuations of V n (t, r) are of order n 1 / 4 and similar to the current fluc- 
tuations of classical walks from the previous section. Recall the definitions 
of Ti and T 2 from (2.10)-(2.11) and abbreviate 

r((*,g),(t,r)) =IMiT 1 ((a,q),{t,r)) + a%T 2 ((s,q),(t,r)). (3.7) 
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Let (V,Z) = (V(t,r),Z(t) : t £ R+,r g R) be the process whose joint 
distribution is defined as follows: 

(i) Marginally, Z(-) = <J2W(-) for a standard Brownian motion W(-). 

(ii) Conditionally on the path Z(-) £ C(R+, R), V is the mean zero Gaus- 
sian process indexed by R + x R with covariance 

E[V(s, q)V(t, r) \ Z{-)] = T((s, q + Z(s)), (i, r + Z(t))) (3.8) 

for (s,q), (t,r) £ R + x R. An equivalent way to say this is to first take 
independent (V°, Z) with Z as above and V° = {V°(t, r) : (t, r) £ R + x R} 
the mean zero Gaussian process with covariance T((s, q), (t, r)J from (3.7), 
and then define V(t, r) = V°(t, r + Z{t)). 

Theorem 3.3. Under the averaged probability P, as n — > oo, the finite- 
dimensional distributions of the joint process {(n~ 1 / 4: V n (t, r), n -1 / 2 Z nt (uj)) : 
t £ R + , r £ R} converge to those of the process (V, Z). 

Thus up to a random shift of the spatial argument we see the same limit 
process as for classical walks: the process V(t, r) = V(t, r — Z{t)) is a mean 
zero Gaussian process with covariance E[V(s,q)V(t,r)} = T[(s,q), (t,r)) 
from (3.7). 

As for classical walks, let us look at the stationary case. The invariant 
distribution is now valid under a fixed ui: the {r] x (0)} are independent and 

(oo i \ 
i=l j=l 

In this case, £ , w ?7o(0) =Var aj 770(0) = /io/(w). By Assumption 1 Ep{p 2 , +£ ) < 
1 for some e > 0, and from that it can be shown that E p[f (uj) 2+e ] < 00. 
Therefore Assumption 2 holds. One can also check that, as for classical 
walks in (2.3), in this stationary situation the flux is linear: H(jj,o) = /JrjVp. 

Recall from Corollary 2.1 that for classical random walks the limit process 
(with fixed space variable r) in the case po = a 1 is fractional Brownian 
motion £ with covariance 

mm)} = ^=(V~ S + vt- vw^\)- 

For RWRE, /jq = o~ 2 implies that 



E[V(s,0)V(t 7 0)\Z(-)} 

= Mo [* CTfs (-^(s)) + * a f t (Z(t)) - * CTf|s _ t| (Z(f) - Z(s))]. 
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Since the right hand side of (3.10) is a non-constant random variable, the 
marginal distribution of V(t, 0) is non-Gaussian. Taking expectations of 
(3.10) with respect to Z(-) gives that 



nV( S ,0)V(t,0)} = ai (Vs + Vt- vU^ff). (3.11) 



Thus we get this conclusion: if fiQ = crfi for RWRE then the limit process 
V (•,()) has the same covariance as fractional Brownian motion, but it is 
not a Gaussian process. 

As the reader may have surmised, we can remove the random shift Z 
from the limit process V by introducing the environment-dependent shift 
in the current process itself. We state this result too. For (t, r) € R+ x M 
define 



(3.12) 



»7m(0) 

m>0 k=l 

r; m (0) 

- X/ zJ 1 i X nt' k > ntv P - Z nt {uj) + ryfn} 

m<0 fe=l 

and its centered version 

V^(t,r) =Y^(t,r) - E u Y^(t,r). 

The process Vn 9 ^ has the same limit as classical random walks. Let V° = 
{V°(t, r) : (t, r) € R + x K} be the mean zero Gaussian process with covari- 
ance r((s, q), (t, r)) from (3.7). 

Theorem 3.4. Under the averaged probability P, as n — > oo, the finite- 
dimensional distributions of the joint process {(n -1 / 4 V^(t, r), n -1 / 2 Z nt (uj)) 
t £ l£+,r € M} converge to those of the process (V°,Z) where V° and Z 
are independent. 



3.2 Sketch of the proof for the quenched mean of the 
current 

The basic thrust of the proof of Theorem 3.3 for the limit of the centered 
current is similar to the one outlined in Section 2.2 for classical random 
walks. The differences lie in the technical details needed to handle the 
random environment. So we omit further discussion of that theorem and 
of the related Theorem 3.4. In this section we explain the main ideas 
behind the proof of Theorem 3.2 for the quenched mean of the current. To 
put aside inessential detail we drop the uniformity, fix (t,r), and sketch 
informally the argument for the following simplified statement: 
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Proposition 3.1. For any e > 0, 

lim ?(w : \E u Y n (t,r) + fj, ry/n + ^ Z nt (uj)\ > £\/n\ = 0. (3.13) 

This proceeds via a sequence of estimations. Abbreviate the centered 
quenched mean by 

W n = E u Y n (t, r) + norVn 

= E u [rio(!n)]P u (X% > ntv P + ry/n) 



m<0 



E u [r] {m)]P u {XZ < ntv P + ry/n) + ^V™. 



m>0 

For a suitable sequence a(n) /* oo define 
W n = £ EMm))Kft ( - f! ^ 

m— — \_a{n)y/n\ +1 



Z nt {9 m u;)-m 



El Z n A9 m uj) — m 
EMm))^ a 2 t ( " n jd + r J + ^ry/n. 



m— 1 



The quenched CLT (part 2 of Theorem 3.1) implies that the difference 



sup 

x6B 



X n - nv P + Z n (u)) 



< 



vanishes P-a.s. as n — ¥ oo. Thus it is possible to choose a(n) /* oo slowly 
enough so that 

-4= W n -W n 



> e = 0. 



lim P 



The next step is to remove the shifts from Z nt by defining 



(3.14) 



m— — la(n)*/nj +1 

[a(n)- N /nJ 



Z nt {oS) - rn 



Z nt (u) - m 



: h r + nory/n. 



The estimation 



lim P 



W n -W n 



> £ = 



(3.15) 



follows from representation (3.1) of Z n t as a sum of ergodic terms whose 
behavior is well understood. 
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Subsequently the ergodicity of the environment allows us to average the 
quenched means E u (r]o(m)) and replace W n with 

ut ^ ( Z nt (ui)-m \ 

w n = y, { — r ) 

m=— la(n)^/n] +1 

|_a(ra)v/nj . . 

El Z n t(ui) — m \ r- 
( + r I + /i rV^- 

m— 1 \v / 

After this, replace the sums with integrals to approximate n _1 / 2 W n with 

Mor - J ^ajt y + r - - $ ff 2 f ^ — r - a; J da:. 

A calculus exercise shows that 




<JV (z - x) - $ Q 2 (-2 - x) dx = z + * Q 2 (A + z) - * Q 2 (A - z), 



where ^ a 2(x) is again the function defined in (2.9). Therefore, 

Jo ^{— +r ~ X )-^{-—- r - X ) dX 

V" V v n / 

- (a(») - — rj . 

The last two terms vanish as n — > oo because ^ CT 2 t (oo) = and 'nr 1 ^ 2 Z nt 
is tight by part 3 of Theorem 3.1. 

Collecting the steps (and taking the omitted details on faith) leads to 
the estimate 

W n = — Mo-^n* + o(\/n) in probability 
as was claimed in (3.13). 
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4 Random average process 
4.1 Model and results 

The state of the random average process (RAP) is a height function a : 
Z — > M where the value <r(i) can be thought of as the height of an interface 
above site i. The state evolves in discrete time according to the following 
rule. At each time point s = 1, 2, 3, . . . and at each site k £ Z, a random 
probability vector Wfc )S = (uik,s(j) ■ —R < i < -R) of length 2i? + 1 is drawn. 
Given the state cr s -i = (er s _i(i) : i £ Z) at time s— 1, at time s the height 
value at site fc is updated to 

a s (k)= ^k, s {j)<y s -i(k+ j). (4.1) 

j:\j\<R 

This update is performed independently at each site k to form the state 
a s = (cr s (k) : k £ Z) at time s. The weight vectors {wfe jS }fe 6 z jS6 N are 
i.i.d. across space-time points (k, s). This system was originally studied by 
Ferrari and Fontes [FF98]. 
Let 

p(0,j)=E[w Ol0 (j)] 
denote the averaged weights with mean and variance 

V = J2xp(x) and of =^2(x - V) 2 p(0,x). (4.2) 

a; xGZ 

Let b = -V. 

Make two nondegeneracy assumptions on the distribution of the weight 
vectors. 

(i) There is no integer h > 1 such that, for some x £ Z, 

^p(0,x + fc/i) = 1. 

feez 

This is also expressed by saying that the span of the random walk with 
jump probabilities p{0,j) is 1 [Dur04, page 129]. It follows that the additive 
group generated by {x £ Z : p(0, x) > 0} is all of Z, in other words this 
walk is aperiodic in Spitzer's terminology [Spi76]. 

(ii) Second, we assume that 

P{maxw .o(7) < 1} > 0. (4.3) 

j 

Let a a be a random average process normalized by 0o(O) = and whose 
initial increments {r)i(0) = <Jo(i) — 0o(i — 1) : i £ Z} are i.i.d. such that 

there exists a > such that E[ \Vi(0)\ 2+a } < oo. (4.4) 
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As before, the mean and variance of initial increments are 

Mo = E{ Vi (0)) and a 2 = Var(%(0)). 

The initial increments 770 are independent of the weight vectors {uj k 

Again we study a suitably scaled process of fluctuations in the charac- 
teristic direction: for (t, r) G R+ x R, let 

Y n (t,r) = n-^ 4 {a^ ntj ([rVE\ + [ntb\) - nory/n}. 

In terms of the increment process 

Vi{s) = <y s (i) - a s {i - 1), 

Y n (t,r) is the centered and scaled net flow from right to left across the 
path s 1 y \jy/n J + \ nsb\ , during the time interval < s < t, exactly as for 
independent particles. 

Recall the definitions (2.10) and (2.11) of the functions Ti and IV 

Theorem 4.1. Under the above assumptions the finite- dimensional dis- 
tributions of the process {Y n (t,r) : (t,r) G R + x R} converge weakly as 
n — > 00 to the finite- dimensional distributions of the mean zero Gaussian 
process {Z(t,r) : (t,r) S M + x R} specified by the covariance 

EZ(s, q)Z(t, r) = ^kT^s, q), (t, r)) + a 2 a T 2 ((s, q), (t, r)). (4.5) 



The constant k is determined by the distribution of the random weights 
and will be described precisely later in equation (4.29). 

Invariant distributions for the general RAP are not known. The next 
example may be the only one where explicit invariant distributions are 
available. 

Example 4.2. Fix positive real parameters 9 > a > 0. Let {ojfc. s ( — 1) : 
s G N, k € Z} be i.i.d. Beta(a, 9 — a) random variables with density 

Hu)= - ( y*) u °-i(i- u) «-°-i 

I (a)l (9 — a) 

on (0,1). Set Wfc jS (0) = 1 — cjfe. s ( — 1). Thus the weights are supported 
on { — 1,0}. A family of invariant distributions for the increment process 
j](s) = (i]k(s) ■ k G Z) is obtained by letting the variables {i]k ■ k G Z} be 
i.i.d. Gamma(#, A) distributed with common density 

fix) - ^Ae-^Az) 9 - 1 (4.6) 

on R + . This family of invariant distributions is parametrized by < A < 
00. Under this distribution E x {q k ] = 9/X and Var A [?; fc ] = 9/X 2 . In this 
situation we find again the fractional Brownian motion limit: 

EZ(s, 0)Z(t, 0) = ci (Vs + Vi - y/\t - s\ ) . (4.7) 

for a certain constant c\. 
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4.2 Steps of the proof 

1. Representation in terms of space-time RWRE 

Let oj = (w/c.s ■ s G N, k £ Z) represent the i.i.d. random weight vectors 
that determine the dynamics, coming from a probability space (f2,6,P). 
Given w and a space-time point (i, r), let {X l s T : s € Z + } denote a random 
walk on Z that starts at Xq T = i and whose transition probabilities are 
given by 

= y | = x) = u x , T - s (y - x). (4.8) 

P w is the path measure of the walk X] ,T , with expectation denoted by E" . 
Comparison of (4.1) and (4.8) gives 

a s (i) = J2P U1 (Xi' 8 =j\Xo' S =i)° S -i(j) =E»[a s . 1 (Xi> s )]. (4.9) 

3 

Iteration and the Markov property of the walks Xl' s then lead to 

a a (i)=E"[a (Xi> s )}. (4.10) 

Note that the initial height function ao is a constant under the expectation 
E u . 

Let us add another coordinate to keep track of time and write X s ' = 
(XI' t ,t — s) for s > 0. Then X^f is a random walk on the planar lattice 
Z 2 that always moves down one step in the e2-direction, and if its current 
position is (x, n), then its next position is (x + y,n — 1) with probability 
u x ,n{y — x). We could call this a backward random walk in a (space-time, 
or dynamical) random environment. 

The opening step of the proof is to use the random walk representation to 
rewrite the random variable Y n (t, r) in a manner that allows us to separate 
the effects of the random initial conditions from the effects of the random 
weights. Abbreviate 

y(n) = [ntb\ + [ry/n\. 
and recall = Eiy(0) and cto(0) = 0. 



a [nti (y(n))=E»[a (Xltt ) i ' lnti )] 



x 



E u 



v(n), |n*J 

L»*J 



1 {X»W'L-*J>0} E ^(°) 





E 



- X l nt ] + 1 



E ^(°)^ ^ *S } ' L " tJ } - E {* > lnti } 



i>0 i<0 

(i H n (t,r) + S„(t,r) 
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where 

H„( ij r)=^(l{i>0}P-{i<X^Lntj } 

- l{^<0}P"{^>X[' I ( ;; J U " tJ }) 

and 

Sn(t, r) = Mo) (l{i > 0}P"{i < xgf L " tJ } 

At this point the terms H n and S n are dependent, but in the course of 
the scaling limit they become independent and furnish the two independent 
pieces that make up the limiting process Z. The limits n _1 ' 4 (i?„(i, ?") — 

r^fn) H(t,r) and n~ 1 / i S n (t,r) S{t,r) are treated separately, and 
then together 

Y n = rT x ' A {H n -ryfii+S n ) A H + S= Z 

with independent terms H and S. This independence comes from the 
independence of the initial height function <to an( l the random environment 
u) that drives the dynamics. The idea is represented in the next lemma. 

Lemma 4.1. Let r\ and uj be independent random variables with values 
in some abstract measurable spaces. Let h n (w) and s n (w,r)) be measur- 
able functions of (to, n). Let E u (-) denote conditional expectation, given ui. 
Assume the existence of random variables h and s such that 

(i) h n (u) h; 

(li) for all 9 € K, E u [e i6an ] -> E(e l0s ) in probability osn^oo. 
Then h n + s n — > h + s, where h and s are independent. 
Proof. Let 9, A € M. Then 

\E(E u [e lXh " +t9s "}) - E[e lXh ] E[e i9s ]\ 

< \E [e lXh " (E u e ie " n - Ee l9s )} \ + \ (Ee lXh " - Ee lXh ) Ee l9s \ 

< \E [e lXh " (E u e i9sn - Ee l9s )] | + \Ee lXh " - Ee lXh \ . 

By assumption (i), the second term above goes to 0. By assumption (ii), the 
integrand in the first term goes to in probability. Therefore by bounded 
convergence the first term goes to as n — > 00. □ 
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We discuss the term S n briefly and reserve most of our attention to H n . 
Two limits combine to give the result. The idea is to apply the Lindebcrg- 
Feller theorem to S n (t,r) under a fixed ui. Then the w-dependent coeffi- 
cients provide no fluctuations but instead converge to Brownian probabil- 
ities due to a quenched central limit theorem for the space-time RWRE. 
Here is an informal presentation where we first imagine that the coefficients 
can be replaced by deterministic quantities: 



S n (t,r) = ^2(vx{0) - mo 
x ( l{x > 0}P 



T^L™ tb J + L''VnJ , L™*J _ 

[nt\ V 



X 

> — - r 



. jr[ntb\ + lr^/n\ , \ nt\ _ r- 

l{x < 0}P"<j = — <^-r 

\/n Jn 



EM ) - W)) > 0}P{B ajt > JL - r} 
-l{x<0}p{B ajt <^-r} 



.re 



Now apply the Lindeberg- Feller theorem to the remaining sum of indepen- 
dent initial occupation variables, and the limiting covariance comes as: 



- 1 / 2 J2^E[S n (t t ,r l )S n (t J ,r J )} 

"Si: 9 *»"" 2 [i; p {^.. > ^-'■•} p { B 'i; > ^- r '} 

x < 
/>oo 

4 Y\ / P^t, >x- n}P{B a 2 t] >x- r 3 }dx 

+ [ P{B a 2 U <x- n}P{B a 2 t . <x- rj } dx 

J — OO 

^o 2 E^ 3 r2(((,r,) ! fe,r J )). 
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Turning this argument rigorous gives the limit in P-probability 

lim E w \e i 'Ef< 0kSn ( tk > rk ^] = E\e i ^ k 9 * s ( t '» r 0l 

n— >oo ^ 

= cxp{-ia 2 ^^ j r 2 ((t l ,r l ), (tj,^))}. 



(4.11) 

i,3 



In particular, in the limit the fluctuations of S come from the initial 
occupation variables r\i (0) and hence are independent of the weights w that 
determine H n . 

2. Quenched mean of the backward space-time RWRE 

The remaining piece of the fluctuations comes from 

H n (t,r) = ^-V^^L^J+L^J, L»*J _ (4 12) 

Theorem 4.3. In the sense of convergence of finite- dimensional distri- 
butions, H n — > H where H(t, r) is the mean zero Gaussian process with 
covariance 

E-ff (s, q)H(t, r) = kTi((s, q), (t, r)) 

« r ?(t+s) l r l , ,1, ( 4 - 13 ) 

dv. 



The constant k is defined below in (4.29). 

By comparing covariances (Exercise 2.1) one checks that H can also be 
defined by 



H(t,r) = ^ y al{t _ s) {r-z)dW{s,z). (4.14) 

J J [0,t] xHL 

Formula (4.14) implies that process {H(t,r)} is a weak solution for this 
initial value problem of the stochastic heat equation: 

H t = ^H rr + y/HW, #(0,r) = 0. (4.15) 



Proof of Theorem 4.3 happens in two steps: first for multiple space points 
at a fixed time, and then across time points. 

Step 1. Martingale increments for fixed time. 

Let us abbreviate X* = x 5 L " tbJ + L ' v " :j - L "* J . Note that 

E ( x *nt\) = l ntb \ + l r Vn~\ + l nt \ v = rVn + 0(l) 
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so H(t,r) in (4.12) is essentially centered and we can pretend that it is 
exactly centered. Let 



g{u) = E"{X Q l '°)-V 



be the centered local drift. Recall the space-time walk X m = (X^ s ,s— m). 
By the Markov property of the walk 



E^(Xr) -x-nV = Y J E" [Xft x - Xl & - V] 

k= ° (4.16) 

n-l n-1 v ' 

= Y^E" [E T <*l"> u (X°<° - V)] = £ E^giT^). 

k=0 k=0 



{T Xjm u})y tS = u x+v ^ m+s is the space-time shift of environments. The g- 
terms above are martingale increments under the distribution P of the 
environments, relative to the filtration defined by levels of environments: 
writing (!>,„,„ = {u> XtS : x € Z, m < s < n}, and with fixed (x, m) and time 
n = 0,1,2..., 



E [£ w g(T^,™w) | U m -n+l,m] 

= P^ixT = (y,m- n)} J g{T y „uj) P(dw m _„) = 0. 



The point above is that the probability P"{X^' m = (y,m — n)} is deter- 
mined by u) m - n+ i. m . 

It turns out that we can apply a martingale central limit theorem to 
conclude that, for a fixed t, a vector 



(ff„(i,ri), H n (t,r 2 ), ■ ■ .,H n (t,r N )) 



becomes a Gaussian vector in the n — > oo limit. Let us take this for 
granted, and compute the covariance of the limit. This leads us to study 
an auxiliary Markov chain which has been useful for space-time (and more 
general ballistic) RWRE. 

Given points (t,q) and (t,r), abbreviate X^ } = x jr ntbi + [q ^ i ^ nti and 
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j£-(2) _ j£[ntb\ + \rs/ri\, [_nt\ 



E[H n (t,q)H n (t,r)] 



^«lj-L^J)^(^j-L^J) 



(2) 



"1/2 E 



[ntj - 1 



[ratj-l 



k=0 



= n-^ 2 J2J2 E [ pU1 ( X ? } = *) p "(4 2) = y)g(T mAnti - j w)g(T yAnti _ kU ) 

j.k x,y 
\nt\-l 

fe=0 

The last step uses the independence of the environment. We denote the 
variance of the drift by <j 2 d = E(g 2 ). Let now Y}. = Xj? — x9~' be the 
difference of two independent walks in a common environment. Yj. is a 
Markov chain on Z with transition probability 



q{x,y) 



^E[w o ,o(z)w o ,o(0, z + y)} x = 

zez 

^p(0,z)p(0,z + 2/-a;) x + 0. 



y,! can be thought of as a symmetric random walk on Z whose transi- 
tion has been perturbed at the origin. The corresponding homogeneous, 
unperturbed transition probabilities are 

q(x, y) = q(0, y-x) = ^2p(0, z)p(0, z + y-x) (x, y e Z). 



Continuing from above, with x n = lr^/n\ — [q^/nl, 
2 L»tJ-i 



E 



[H n {t,q)H n {t,r)]=^ V q k {x n , 0) = ^G [nti -i(x n , 0). (4.17) 



If Yfe were a symmetric random walk, we would know this limit exactly 
from the local central limit theorem: 

Lemma 4.2. For a mean 0, span 1 random walk S„ on Z with finite 
variance a 1 , bgM and points a n £ Z such that \a n — a,y/n\ = 0(1), 



[ntj-l 



lim "4? II P{S k = a n ) = \ [ — =^cxp{-^-| 
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Proof. By the local CLT [Dur04, Section 2.5] 



lim sup a/to P{S m = x) 1 = exp{- — - | 



0. 



v27rm<7^ 

Use this in a Riemann sum argument (details as exercise). □ 
For the homogeneous q-walk this lemma gives (using symmetry) 

1 L'^ 1 1 f 2 ^ 1 r (r-o) 2 ^ 

lim V" qk( x )= -=cxp\- { —--^- \ dv. (4.18) 

Now the task is to relate the transitions q and 5. For this purpose we 
introduce one more player: the potential kernel of the symmetric g-walk, 
defined by 

a(x) = lim [G n (0,0)-G„(x,0)] 

n— >-oo 

= lim ^g- fe (0,0)-^ g - fc (x,0) . 

k fe=0 k=0 ' 

The potential kernel satisfies a(0) = 0, the equations 



e) =^2q(x,y)a(y) for x / 0, and £ g(0, y)a(y) = 1, (4.20) 
ye 

and the limit 



a(x) 



1™..^ = ^. (4.21) 



a(x) 1 

a;->±oo |x| 2(7^ 

(For existence of a and its properties, see [Spi76, Sections 28-29].) 



Example 4.4. If for some k <G Z, p(0, fc)+p(0, fc+1) = 1, so that q(0, x) = 
for x ^ {-1, 0, 1}, then a(x) = |a:[/(2cr|). 

Define the constant 

/3 = 5>(0,x)a(x). (4.22) 

This constant accounts for the difference in the limits of the Green's func- 
tions for transitions q and q. 

Lemma 4.3. Let ifM and x n £ Z be such that x n — n 1 ' 2 x stays bounded. 
Then 

lim n- 1 / 2 G„(x„,0) = — ^ / ' -rt= expj-^-) dv. (4.23) 
n^oo 2/30- j J \/2irv L 2w J 
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Proof. Given limit (4.18), it suffices to prove 



sup 



4=G„(z,0) - -^G„0,0) 



0. 



First we prove the case z = 0. 
By (4.18), 



We need to show 



lim n~ 1/2 G n (0,0) 



lim n" 1/2 G„(0,0) 



(4.24) 

(4.25) 
(4.26) 



Using (4.20), a(0) = 0, and q(x,y) = q(x,y) for x ^ 0, 

Y, 9 m {0, x)a(x) = Y 9 m (0, x)a(x) = £ g m (0, x)q(x, y)a(y) 

= ]T g m (0, 

= Y g m+1 (0, y)a(y) - g m (o, o) £ g(0, y)a( y ). 

t/GZ yGZ 

Constant /3 appears in the last term. Sum over m = 0, 1, . . . , n — 1 to get 
(1 + q(0. 0) + ■ ■ ■ + q n -\0, 0))/3 = ^ g"(0, a;)o(x) 

and write this in the form 

n- 1/2 /3G n _!(0,0) = n- 1 / 2 £; [a(^n)]. 

Recall that y n = X n — X n where X n and X n are two independent walks 
in the same environment. By the quenched CLT for space-time RWRE, 
n~ r ^Y n A/"(0, 2tr 2 ). Marginally X n and X n are i.i.d. walks with 

bounded steps, hence there is enough uniform integrability to conclude 
that 

n- l ' 2 E a \Y n \^2^al/^. 
By (4.21) and some estimation (exercise), 

1 



n-^ 2 E [a{Y n )] 



This proves (4.26) and thereby limit (4.24) for 2 = 0. 

To get the full statement in (4.24), for k > 1 and z ^ let 

f k (z,0) = !{ 2 #o} Y q{z,zi)q{z u z 2 )---q{z k -i^) 

Z1^0,...,Zfe_i#0 



(4.27) 



Current fluctuations 



31 



denote the probability that the first visit to the origin occurs at time k. 
This quantity is the same for both q and q because these processes do not 
differ until the origin is visited. Choose no so that 



|/3G m (0,0) — G m (0, 0)| < £\fm for m > n a . 



Then 



sup 



-£=G n (z,0) - -^G„(z,0) 
\ n Jn 



1 ™ 

<su P — V/ fc (z,0)|/3G„^(0,0) - G„_ fc (0,0) 

< sup JL y / fe (z ,o)v / ^fc + ^><e+^ 
V"- ^ V" v« 



Letting n — > oo completes the proof. □ 
Combining (4.17) and (4.23) gives 



lim E[H n (t,q)H n (t,r)] = f ' -^exp{ 

n->oo 2^crf Jo \/2irv I- 



2/3<7? 7 v 7 ^ 1 2w J (4.28) 
= Kri((t,«),(t,r)), 
where wc defined a the new constant 



While we have not furnished all the details, let us consider proved that 
for a fixed t, the finite-dimensional distributions of H(t,r) converge to the 
Gaussian process H(t,r) with covariance K,Ti((t,q), (t,r)). 

Step 2. Markov property for time steps. 



This step is overly technical and so we only give a sketch of the idea 
ehind it. 
[ns\ gives: 



behind it. Stopping and restarting the walk ^L ntb J + L'^J, L n *J ^ i eve ] 



' l_ n *J — L ns J J ' / V |_ ns J 

+ E"{xW+£f*'M) - \nsb\ \ryfr\. 
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Change summation index to u = x/^/ri. Then we have approximately the 
identity 



where H*(t — s, r) is the same as H n (t — s, r) but with origin shifted (ap- 
proximately) to (nsb,ns). On the right-hand side, the processes H n (s, •) 
and H*(t — s, •) are independent of each other because they depend on 
disjoint levels of environments: H n (s, •) uses u>i i ns j and H*(t — s, •) uses 
^InsJ+i, \nt\ ■ As n — >■ oo the probability coefficients of the sum converge to 
deterministic Gaussian probabilities by the quenched CLT for the RWRE. 
By the result for fixed t, the right-hand side above converges in distribution. 
Taking the limits and supplying all the technicalities leads to the equation 



where on the right, the processes H(s, •) and H*(t— s, •) arc independent. 
From this equation one can verify that the finite-dimensional distributions 
of the process H(t,r) are Gaussian with covariance nTi((s, q), (t,r)) as 
stated in Theorem 4.3. 

This concludes the presentation of the random average process limit. 

References 

The fluctuation results for RAP presented here are from [BRAS06]. In 
addition to [FF98], RAP was later studied also in [FMV03] . The quenched 
CLT for space-time RWRE has been proved several times with progressively 
better assumptions, sec [RAS05]. 

5 Asymmetric simple exclusion process 

The asymmetric simple exclusion process (ASEP) is a Markov process that 
describes the motion of particles on the one-dimensional integer lattice Z. 
Each particle executes a continuous-time nearest-neighbor random walk on 
Z with jump rate p to the right and q to the left. Particles interact through 
the exclusion rule which means that at most one particle is allowed at each 
site. Any attempt to jump onto an already occupied site is prevented from 
happening. The asymmetric case is p ^ q. We assume < q < p < 1 and 
p + q = 1. 



H n (t,r) 




+ K(t 




Current fluctuations 



33 



For this process we do not derive precise distributional limits for the 
current, but only bounds that reveal the order of magnitude of the fluctu- 
ations. In contrast with the earlier results for linear flux, the magnitude of 
current fluctuations is now rj 1 / 3 . 

The proofs of these bounds are based on couplings, and make heavy use 
of the notion of second class particle. 

5.1 Basic properties 

We run quickly through the fundamentals of (p, g)-ASEP. 

Definition and graphical construction. The state of the system at 
time t is a configuration n(t) = (T]i(t)) ie z G {0, 1} Z of zeroes and ones. The 
value r]i(t) = 1 means that site i is occupied by a particle at time t, while 
the value iji(t) = means that site i is vacant at time t. 

The motion of the particles is controlled by independent Poisson pro- 
cesses (Poisson clocks) {N^ i+1 , N^*' 1 : i G Z} on R + . These Poisson 
processes are independent of the (possibly random) initial configuration 
?y(0). Each Poisson clock N 1 ^^ 1 has rate p and each j\H ,_H has rate q. 
If t is a jump time for 7\r>^ l+1 and if (i]i(t— ), 7/4+1 (i— )) = (1,0) then at 
time t the particle from site i moves to site i + 1 and the new values are 
(r)i(t),r]i+i(i)) = (0, 1). Similarly if t is a jump time for TV 4 " 14 " 1 a particle 
is moved from i to i — 1 at time t, provided the configuration at time t— 
permits this move. If the jump prompted by a Poisson clock is not permit- 
ted by the state of the system, this jump attempt is simply ignored and 
the particles resume waiting for the next prompt coming from the Poisson 
clocks. 

This construction of the process is known as the graphical construction 
or the Harris construction. When the initial state is a fixed configuration 
T), P v denotes the distribution of the process. 

We write 77, u>, etc for elements of the state space {0, 1} Z , but also for 
the entire process so that 77-process stands for {t]i(t) : i <E Z, < t < 00}. 
The configuration 5{ is the state that has a single particle at position i but 
otherwise the lattice is vacant. 

Remark 5.1. When infinite particle systems such as ASEP are constructed 
with Poisson clocks, there is an issue of wcll-definedness that needs to be 
resolved. Namely, if we ask whether site x is occupied at time t, we need 
to look backwards in time at all the possible sites from which a particle 
could have moved to x by time t. This might involve an infinite regression: 
perhaps there is a sequence of times t > t\ > t 2 > t 3 > • • • > such that 
Poisson clock Tyz-fc-s-z-fc+i s jg na i e( j a jump attempt at time tk- Such a 
sequence of jumps could in principle bring a particle to x "from infinity." 

However, it is easily seen that this happens only with probability zero. 
For any fixed T < 00 there is a positive probability that both 7V i_M+1 and 
are em pty m jo^ T}. Consequently almost surely there are infinitely 
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many edges (i, i + 1) across which no jump attempts are made before time 
T. This way the construction can actually be performed for finite portions 
of the lattice at a time. 

Similar issue arises with the possibility of simultaneous conflicting jump 
commands. By excluding a zero-probability set of realizations of {N 1 ^^ 1 } 
we can assume that there are no simultaneous jump attempts. 

Invariant distributions. A basic fact is that i.i.d. Bernoulli distri- 
butions {^ p } p e[o,i] ar c extremal invariant distributions for ASEP. For each 
density value p £ [0, 1], v p is the probability measure on {0, 1} Z under which 
the occupation variables {rji} are i.i.d. with common mean fr]idi> p = p. 
When the process ry is stationary with time-marginal v p , we write P p for 
the probability distribution of the entire process. The stationary dcnsity-p 
process means the ASEP 77 that is stationary in time and has marginal 
distribution r](t) ~ v p for all t G M + . 

Remark 5.2. A note about how one would check the invariance. In general, 
the infinitesimal generator L of a Markov process is an operator defined as 
the derivative of the semigroup: 

Above E n denotes expectation under P n , the distribution of the process 
when the initial state is rj. The generator of ASEP is 



(5.2) 



it: 



that acts on cylinder functions ip on the state space {0, 1} Z and rf^ ~ 
77 — Si + dj is the configuration that results from moving one particle from 
site i to j. Equation (5.1) can be derived from the graphical construction 
with some estimation. 

Invariance of a probability distribution can be checked by a generator 
computation. For ASEP it is enough to check that 

Lip dp = (5.3) 

for cylinder functions <p to conclude that p is invariant. This can be used 
to check that Bernoulli measures v p are invariant for ASEP. 

Basic coupling and second class particles. The basic coupling of 
two exclusion processes rj and oj means that they obey a common set of 
Poisson clocks {A^~" +1 , N 1 ^ 1 ^ 1 }. Suppose the two processes 77 and r} + 
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satisfy rj + (0) = 77(0) + <Wo) at time zero, for some position Q(0) £ Z. This 
means that 77+ (0) = r)i(0) for all i ^ Q(0), ?7q( )(0) = 1 and ?7q( )(0) = 
0. Then throughout the evolution in the basic coupling there is a single 
discrepancy between rj(t) and ij + (t) at some position Q(t): r/ + (t) = rj(t) + 
(5g( t ). From the perspective of r](t), Q(t) is called a second class particle. 
By the same token, from the perspective of f] + (t), Q{t) is a second class 
antiparticle. In particular, we shall call the pair (ij, Q) a (p, q)-ASEP with 
a second class particle. 

We write a boldface P for the probability measure when more than 
one process are coupled together. In particular, P p represents the situ- 
ation where the initial occupation variables ^i(O) = rj^~(0) are i.i.d. mean-p 
Bernoulli for i 7^ 0, and the second class particle Q starts at Q(0) = 0. 

More generally, if two processes rj and co are in basic coupling and w(0) > 
77(0) (by which we mean coordinatewise ordering Wj(0) > ?7 ? :(0) for all i) 
then the ordering uj{t) > rj(t) holds for all < t < 00. The effect of 
the basic coupling is to give priority to the 77 particles over the lo — n 
particles. Consequently we can think of the cj-process as consisting of 
first class particles (the -q particles) and second class particles (the lo — r\ 
particles) . 

Current. For x £ Z and t > 0, J x (t) stands for the net left-to-right 
particle current across the straight-line space-time path from (1/2,0) to 
(x + 1/2, t). More precisely, J x (t) = J x (t) + — J x (t)~ where J x {t) + is the 
number of particles that lie in (—00, 0] at time but lie in [x+ 1, 00) at time 
t, while J x (t)~ is the number of particles that lie in [1, 00) at time and in 
(—00, x] at time t. When more than one process (w, 77, etc) is considered in 
a coupling, the currents of the processes are denoted by J x (t), J x (t), etc. 

5.2 Results 

The average net rate at which particles in the stationary (p, g)-ASEP at 
density p move across a fixed edge (i, i + 1) is the flux 

H(p) = E'Vo(t)} = (p - q)p(l - P). (5.4) 

This formula follows from the fact that this process M{t) is a mean zero 
martingale: 

M(t) = J Q (t) - f (pl{ m (s) = 1, mis) = 0} 

Jo v (5.5) 

- gl{?7i(s) = 1, ??o(s) = 0}) ds. 

For the more general currents 



E p [J x {t)] = tH{p) - xp (i£Z,(>0) 



(5.6) 



36 



Timo Seppalainen 



as can be seen by noting that particles that crossed the edge (0,1) either 
also crossed (x, x + 1) and contributed to J x (t) or did not. 
The characteristic speed at density p is 

V = H'(p) = (p-q)(l-2p). (5.7) 

The derivation of the fluctuation bounds for the current rests on several 
key identities which we collect in the next theorem. 

Theorem 5.3. Let the second class particle start at the origin: Q(0) = 0. 
For any density < p < 1, z e Z and t > we have these formulas. 

Var"[J 2 (<)] = £|j - z\ C O v'[Tfe(t) > tto(0)], (5.8) 

jez 

Cov'te(t), tjo(0)] = p(l - P )P p {Q(t) = j}, (5.9) 

and 

E p [Q(i)] = V p t. (5.10) 



Formulas (5.8) and (5.9) combine to give 

Var"[J,(t)] = p(l - p)W\Q{t) - z\. (5.11) 

In particular, for the current across the characteristic, 

V^[J\y Pt] {t)]=p{l-p)W\Q{t)-[Vn\\. (5.12) 

Thus to get variance bounds on the current, we derive moment bounds on 
the second class particle. 

We now state the main result, the moment bounds on the second class 
particle. It is of interest to see how the bounds depend on the bias 9 = p — q 
so we include that in the estimates. 

Theorem 5.4. There exist constants < Cq,C < oo such that, for all 
< 9 < 1/2, < p < 1, 1 < m < 3, and t > c 9~ 4 , 

K2m/3 t 2m/3 < £ p r ,q,q _ yp^m j < ^^^m/S^m/S^ (5.13) 

C 3 — m 

For the upper bound the constants are fixed for all values of the parameters. 
For the lower bound both constants Cq,C depend on the density p. 

As a corollary for m = 1, we obtain the bounds for the variance of the 
current seen by an observer traveling at the characteristic speed V p : for 
* > c (p)9-\ 

Ci(p)fl 1/3 £ 2/3 < Var p [J LyPtJ (t)] < C 2 1/3 i 2/3 . (5.14) 
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It follows from the variance bound (5.14) that for v ^ V p a Gaussian 
limit in the central limit scale holds: 



J [tv] (t) - EP(J [tv] {t)) 

t l/2 



V 



X (5-15) 



for a centered normal random variable x- To observe this, take the case 
v > V p . Let J* be the current across the straight-line space-time path 
from ((v — V p )t, 0) to (vt, t). This current has variance of order i 2 / 3 . Then 
use 

(v-V)t 

r = j [tv] (t)+ ? ^°)- 

A distributional limit exists for the current for the case of the stationary 
totally asymmetric simple exclusion process (TASEP). We state the result 
here. In TASEP particles march only to the right (say), and so p = 1 and 
q = 0. 

Theorem 5.5. [FS06] In stationary TASEP, the following distributional 
convergence holds: 

lim pP { Avt\ (*) ~ P 2 t 1 Fn(x) (5 16) 

A™ \ f?/3(l - p)2/3il/3 " X j r °^ X) ^- W> 



The distribution function Fo above is defined in [FS06] as Fo(x) = 
(d/dx)(FG\jE{x)g(x,0)) where Fgue is the Tracy- Widom GUE distribu- 
tion and g a certain scaling function. 

Theorem 5.5 will not be discussed further, and we turn to proofs of 
Theorem 5.3 and Theorem 5.4. In the next section we give partial proofs 
of the identities in Theorem 5.3. Section 5.4 describes a coupling that we 
use to control second class particles, and a random walk bound that comes 
in handy. The last two sections of this chapter prove the upper and lower 
bounds of Theorem 5.4. 

5.3 Proofs for the identities 

Let uj be a stationary exclusion process with i.i.d. Bcrnoulli(p) distributed 
occupations {uji(t)} at any fixed time t. 

Proof of equation (5.8). This is partly a hand- waiving proof. What is miss- 
ing is justification for certain limits. 

To approximate the infinite system with finite systems, for each N € N 
let process u> N have initial configuration 

w i Ar (0)=w 4 (0)l { _^< i < JV} . (5.17) 
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We assume that all these processes are coupled through common Poisson 
clocks. Let J^(t) denote the current in process uj n . 

Let z(0) = 0, z(t) = z, and introduce the counting variables 

= £ w n(*), I N W= E "«(*)■ (5.18) 

n>z(t) n<z(t) 

Then the current can be expressed as 

j?(t) = r?(t)-i»(o) = i»(o)-i»(t), 

and its variance as 

Var Jf (t) = Cov(Z* (t) - if (0), I N (0) - I N (t)) 

= Cov(lf (t),l5(0)) + Cov(Zf (0), I N (t)) 

~ Cov(lf (0), 15 (0)) - Cov(/~ (t), I N (t)). 

Independence of initial occupation variables gives 

Cov(/f(0),/ w (0))=0 

and the identity above simplifies to 

Var Jf (t) - Cov(/f(t),/ w (0)) + Cov(lf (0),/ N (t)) 

-Cov(Zf(i),Z5(t)) 

= E Oov[^(t),«f(0)] (5-19) 

+ E Cov[u£(t),ujZ(0)}-Cov(I»(t),I?(t)). 

k<z 1 m>0 

In the N — > oo limit variables u>f{t) converge (a.s. and in L 2 ) to the 
i.i.d. occupation variables uii(t) of the stationary process. It follows from 
the graphical construction that on a fixed time interval covariances can be 
bounded exponentially, uniformly over N: for a fixed < t < oo, 

| Cov^i),^ (s)] I < Ce- Cl l™- fc l for a € [0,*]. 

Hence in the limit the last covariance in (5.19) vanishes. Furthermore, 
Jf(i) J 2 (i) similarly, so in the limit we get 

Var J z {t)= E Cov[w m (i),cj fc (0)] + ^ Cav[w fc (t),w ro (0)] 

A;<0, m>z k<z, m > 

= 5^|n-«|Cov[w n (t),Wo(0)]. 

nez 

This proves equation (5.8). □ 
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Proof of equation (5.9). This is a straight-forward calculation. 

Cov>j(t),wo(0)] = EP[uj(t)uo(p)] - p 2 
= pEP[ Uj (jt) | w (0) = 1] - pE"[uj{t)\ 
= pfE^it) | w (0) = 1] - pE"[ojj(t)\uj (0) = 1] 

- (l-p)EP[oj j (t)\u (0) = 0] 

= p(l - p)(E"[uj(t)\uo(fi) = 1] - |wo(0) = 0]) 

= p(l p)(W[u+(t)] E>j(t)]) = p(l p)W[u+(t) - uj(t)] 

= p(l-p)PP[Q(t)=j]. □ 

Proof of equation (5.10). Let again oj n be the finite process with initial 
condition (5.17). Let I N = be the number of particles in the 

process u . I is a Binomial(2iV + l,p) random variable. For < p < 1 

d ^ / 2 N+1 



J-E[J?(t)] = f E - P) 2Ar+1 - m ^[^(i)l^ = m] 

= ^F(f = - 2iY 1 +1 ' m )^fW|/ jV ^ m] 

3 P P 

1 





p) 


1 




p(l- 


p) 


1 




p(l- 


p) 



jf(i)(l w -(27V+l)p) 



Covllf(i)-lf (0),/^(0)+F' 



The last equality used Cov[if (0), (0)] = that comes from the i.i.d. 
distribution of initial occupations. The first covariance on line (5.20) write 
directly as 



Cov[^(i),/ N (0)]= J2 CovlwZ®,^' 

k<0, m>z 



The second covariance on line (5.20) write as 

Cov[/f (t)-l£ (0), If (0)] = Cov[/ w (0) -I N (t), If (0)] 

= -Cov[/ w (i),/f(0)] = - 2 C OT [ W ?(t),c4(0)]. 

k<z. m>0 
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Inserting these back on line (5.20) gives 

1 

P(l-P) 



J2 Cov[o£(t),u${0)] 



k<0, m>z 

J2 Cov[ W f(t),^(0)] 

&<z, m>0 

Compared to line (5.19) we have the difference instead of the sum. Integrate 
over the density p and take N — ¥ oo, as was taken in (5.19), to obtain 



E^J z (t)}-E x [J z (t)} 



6(1- 
(E 9 [Q(t)]- z)dQ 



^(j-z)Cov%-(i), Wo (0)] 



d0 



for < A < p < 1. Couplings show the continuity of these expectations: 

E x [J z (t)} -^E"[J z {t)] and E A [Q(i)] -> W[Q(t)] (5.21) 

as A p in (0, 1). Thus the identity above can be differentiated in p. With 
z = and via (5.4) identity (5.10) follows. □ 

5.4 A coupling and a random walk bound 

As observed in (5.10) the mean speed of the second class particle in a 
density- p ASEP is H'(p). Thus by the concavity of if a defect travels on 
average slower in a denser system (recall that we assume p > q throughout). 
However, the basic coupling does not respect this, except in the totally 
asymmetric (p = 1, q = 0) case. To see this, consider two pairs of processes 
(lo + ,lu) and (^ + ,77) such that both pairs have one discrepancy: uj + (t) = 
ui(t) + Sq^j^) and rj + (t) = r)(t) + Sq^u). Assume that ui(t) > r](t). In basic 
coupling the jump from state 
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happens at rate q and results in = i + 1 > i 

In this section we construct a different coupling that combines the basic 
coupling with auxiliary clocks for second class particles. The idea is to 
think of a single "special" second class particle as performing a random 
walk on the process ofui — r] second class particles. This coupling preserves 
the expected ordering of the special second class particles, hence it can be 
regarded as a form of microscopic concavity. 

This theorem summarizes the outcome. 
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Theorem 5.6. Assume given two initial configurations {Ci(0)} and {£i(0)} 
and two not necessarily distinct positions Q^(0) and Q^(0) on Z. Suppose 
the coordinatewise ordering £(0) > £(0) holds, Q^(0) < Q^(0), and d(0) = 
&(0) + 1 for i e {Q c (0),Q ? (0)}. Define the configuration C~(0) = C(0) - 

Then there exists a coupling of processes 



with initial state (£ (0), Q^(0), £(0), Q^(0)) as described in the previous 
paragraph, such that both pairs (C",<3 C ) and are {p,q)-ASEP's with 

a second class particle, and Q^(t) < Q^(t) for all t > 0. 

To begin the construction, put two exclusion processes £ and £ in basic 
coupling, obeying Poisson clocks {N l ^ l±1 }. They are ordered so that 
£ > £. The £ — £ second class particles are labeled in increasing order 
• • • < X m _i(t) < X m (t) < X, n+ i(t) < • • • . We assume there is at least one 
such second class particle, but beyond that we make no assumption about 
their number. Thus there is some finite or infinite subintcrval 7 C Z of 
indices such that the positions of the £ — £ second class particles are given 
by {X m (t) : m G 7}. 

We introduce two dynamically evolving labels a(t), b(t) € 7 in such a 
manner that -Xa(t) (t) is the position of a second class antiparticle in the 
^-process, -X&(t) (t) is the position of a second class particle in the ^-process, 
and the ordering 



is preserved by the dynamics. 

The labels a(t), b{t) are allowed to jump from m to m ± 1 only when 
particle X m ±\ is adjacent to X m . The labels do not take jump commands 
from the Poisson clocks {N l ^ l±1 } that govern (£, £). Instead, the directed 
edges + 1) and (i,i — 1) are given another collection of independent 
Poisson clocks so that the following jump rates are realized. 

(i) If a = b and X a+ i = X a + 1 then 



X a(t) {t) < X b(t) {t) 



(5.22) 



(a, b) jumps to 




(ii) If a = b and X a ^\ = X a — 1 then 



(a, b) jumps to 
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(iii) If a 7^ b then a and 6 jump independently with these rates: 

X a + 1 
X a -l; 

X b + 1 
X b -1. 

Let us emphasize that the pair process (£, £) is still governed by the old 
clocks {N l ~^ l±1 } in the basic coupling. The new clocks on edge {i,i + 1} 
that realize rules (i)-(iii) are not observed except when sites {i,i + 1} are 
both occupied by X-particles and at least one of X a or X b lies in {i, i + 1}. 

First note that if initially a(0) < 6(0) then jumps (i)-(iii) preserve the 
inequality a(t) < b(t) which gives (5.22). (Since the jumps in point (iii) 
happen independently, there cannot be two simultaneous jumps. So it is 
not possible for a and b to cross each other with a (a, b) — s- (a + 1, b — 1) 
move.) 

Define processes C~(*) = C(*) ~ s x aW (t) and £+(t) = + Sx b(t) (t)- In 
other words, to produce £ _ remove particle X a from f , and to produce 
£ + add particle X b to £. The second key point is that, even though these 
new processes are no longer defined by the standard graphical construction, 
distributionwise they are still ASEP's with second class particles. We argue 
this point for (£~,X a ) and leave the argument for (£,Xb) to the reader. 

Lemma 5.1. The pair (£~,X a ) is a (p,q)-ASEP with a second class par- 
ticle. 

Proof. We check that the jump rates for the process ((~,X a ), produced 
by the combined effect of the basic coupling with clocks {7V l-M±1 } and the 
new clocks, are the same jump rates that result from defining an (ASEP, 
second class particle) pair in terms of the graphical construction. 

To have notation for the possible jumps, let denote an empty site, 1 a 
£ - -particle, and 2 particle X a . Consider a fixed pair (i, % + 1) of sites and 
write xy with x, y € {0, 1, 2, } for the contents of sites (i, i + 1) before and 
after the jump. Then here are the possible moves across the edge {i, i+ 1}, 
and the rates that these moves would have in the basic coupling. 



Type 1 


10 - 


-> 01 


with rate p 




01 - 


-> 10 


with rate q 


Type 2 


20 - 


-> 02 


with rate p 




02 - 


-> 20 


with rate q 


Type 3 


12 - 


-> 21 


with rate p 




21 - 


-> 12 


with rate q 



a jumps to 
b jumps to 



a + 1 with rate q if X a+ \ = 

a — 1 with rate p if A Q _i = 

6+1 with rate p if X b +i = 

6 — 1 with rate q if X b -i = 
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Our task is to check that the construction of ((~ , X a ) actually realizes these 
rates. 

Jumps of types 1 and 2 are prompted by the clocks {7V l_M±1 } of the 
graphical construction of (£,£), and hence have the correct rates listed 
above. 

Jumps of type 3 occur in two distinct ways. 

(Type 3.1) First there can be a ^-particle next to X a , and then the rates 
shown above are again realized by the clocks {N l ^ l±1 } because in the basic 
coupling the ^-particles have priority over the X-particles. 

(Type 3.2) The other alternative is that both sites {i,i + 1} are occu- 
pied by X-particles and one of them is X a . The clocks {A^^'* 1 } cannot 
interchange the X-particlcs across the edge {i, i + 1} because in the (£, £)- 
graphical construction these are lower priority ^-particles that do not jump 
on top of each other. The otherwise missing jumps are now supplied by 
the "new" clocks that govern the jumps described in rules (i)-(iii). 

Combining (i)-(iii) we can read that if X a = i + 1 and X a -i = i, then 
a jumps to a — 1 with rate p. This is the first case of type 3 jumps above, 
corresponding to a £~-particle moving from i to i + 1 with rate p, and 
the second class particle X a yielding. On the other hand, if X a = i and 
X a+ i — i + 1 then a jumps to a + 1 with rate q. This is the second case in 
type 3, corresponding to a ^"-particle moving from i + 1 to i with rate q 
and exchanging places with the second class particle X a . 

We have verified that the process ((~,X a ) operates with the correct 
rates. 

To argue from the rates to the correct distribution of the process, we 
can make use of the process (C~,C)- The processes (£~,X a ) and (C~,C) 
determine each other uniquely The virtue of (£~ , £) is that it has a compact 
state space and only nearest-neighbor jumps with bounded rates. Hence 
by the basic theory of semigroups and generators of particle systems as 
developed in [Lig85], given the initial configuration, the distribution of 
the process is uniquely determined by the action of the generator on local 
functions. Thus it suffices to check that individual jumps have the correct 
rates across each edge {i, i + 1}. This is exactly what we did above in the 
language of (C - , X a ). □ 

Similar argument shows that (£, X^) is a (p, g)-ASEP with a second class 
particle. To prove Theorem 5.6 take Q 1 * = X a and = Xi,. This gives 
the coupling whose existence is claimed in the theorem. 

To conclude, let us observe that the four processes (£, £ + , £ ~, () are not 
in basic coupling. For example, the jump from state 
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happens at rate q (second case of rule (i)), while in basic coupling this move 
is impossible. 

As the second point of this section we prove a random walk estimate. 
Let Z(t) be a continuous-time nearest-neighbor random walk on state-space 
SCZ that contains Z_ = {..., -2, -1, 0}. Initially Z(0) = 0. Z attempts 
to jump from x to x + 1 with rate p for x < — 1, and from x to x — 1 with 
rate g for x < 0. Assume p > q = 1 — p and let # = p — g. The rates on 
S* s Z_ need not be specified. 

Whether jumps are permitted or not is determined by a fixed environ- 
ment expressed in terms of {0, l}-valued functions {u(x, t) : x £ S, < t < 
oo}. A jump across edge {x — 1, x} in either direction is permitted at time 
t if u(x,t) = 1, otherwise not. In other words, u(x, i) is the indicator of 
the event that edge {x — 1, x} is open at time t. 

Assumption. Assume that for all x £ S and T < oo, u(x,t) flips 
between and 1 only finitely many times during < t < T. Assume for 
convenience right-continuity: u(x,t+) = u(x,t). 

Lemma 5.2. For all t > and k > 0, 

P{Z{t) < -k} < er 2Sk . 

This bound holds for any fixed environment {u{x, t)} subject to the assump- 
tion above. 

Proof. Let Y(t) be a walk that operates exactly as Z(t) on Z_ but is 
restricted to remain in Z_ by setting the rate of jumping from to 1 to 
zero. Give Y(t) geometric initial distribution 

P{Y(0) = -j} = = (l - ^) (^) ] for j > 0. 

The initial points satisfy Y(0) < Z(0) a.s. Couple the walks through 
Poisson clocks so that the inequality Y(t) < Z(t) is preserved for all time 
< t < oo. 

Without the inhomogeneous environment Y(t) would be a stationary, 
reversible birth and death process. We argue that even with the environ- 
ment the time marginals Y(t) still have distribution tt. This suffices for the 
conclusion, for then 

P{Z(t) < -k] < P{Y(t) < -k} = (q/p) k = cxp(fclogi T f) < e~ 26k . 

To justify the claim about Y(t), consider approximating processes 
m £ N, with the same initial value y( m '(0) = Y(0). Y^ m \t) evolves so 
that the environments {u(x, t)} restrict its motion only on edges {x — 1, x} 
for — m + I < x < 0. In other words, for walk Y^ m ^(t) we set u(x, t) = I 
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for a; < —m and < t < oo. We couple the walks together so that 
Y(t) = Y^ m \t) until the first time one of the walks exits the interval 
{-m+1,. ..,()}. 

Fixing m for a moment, let = sq < s\ < S2 < S3 < . . . be a partition 
of the time axis so that Sj /* 00 and the environments {u{x,t) : — m < 
x < 0} are constant on each interval t G [sj,Sj+i). Then on each time 
interval [sj, s»+i) y' m '(i) is a continuous time Markov chain with time- 
homogeneous jump rates 

Jpu(x + 1, Si), -m < .t < 
c{x,x + 1) = < 

I p, x < —m — 1 



and 

, . \qu(x,Si), -m + l<x<0 

c(x, x — 1) = < 

\q, x < —m. 

One can check that detailed balance tt(x)c(x, x + 1) = ir(x + l)c(x + 1, x) 
holds for all a; < — 1. Thus tt is a reversible measure for walk Y( m )(£) on 
each time interval [si, Si+i), and we conclude that Y^ m '(t) has distribution 
7r for all < t < 00. 

The coupling ensures that Y^ m \t) — > Y"(t) almost surely as m — > 00, and 
consequently also Y(t) has distribution tt for all < t < 00. □ 



5.5 Proof of the upper bound for second class particle 
moments 

Abbreviate 

*(t) = E p \Q{t) - V p t\. (5.23) 

Lemma 5.3. Let B £ (0, 00). Then there exists a numerical constant 
C G (0,oo) and another constant C\{B) G (0,oo) such that, for all densities 
< p < 1, u > 1, < 9 < 1/2, and t > c 1 (B)0- 1 , 

P p {Q(t) > V p t + u} 

t 2 , , , t 



( 1 w t ) + L-) + e~ u2 / c \ B6 2 /H 2 / 3 <u< 20i/3 (5.24) 

\M 4 U 6 J 

e-'^ c , u > 20t/3. 



ce 2 

< J \ W 1 u 



Proof. First we get an easy case out of the way. 
Case 1. u > 5p8t. 
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This comes from an exponential Chebyshev argument. Let Z t be a 
nearest-neighbor random walk with rates p = (1 + 0)/2 to the right and 
q={\-6)/2 to the left. For a g (0, 1], using 



< 1 + a 2 and < a + a 2 , 



we get 



E[e«*] = ex P (-f + t G + e + gt £ 2 6 ) 
< exp(2a 2 < + a9t) . 



(5.25) 



We have the stochastic domination Z t > because no matter what the 
environment next to Q(t) is, Q(i) has a weaker right drift than Z t . Then, 
since V p = 9(1 - 2p) and 2p6t < 2u/5, 



P p {Q(t) > V p t + u} < P{Z t >9t + f u} 

3 
'5 

^ I 



< exp(-faw + 2a 2 t) 



exp(-3u/10) u > 20t/3. 



(5.26) 



In the last inequality above choose a = 1 A Note that 20t/3 > hpOt. 
It remains to consider this range of u: 

Case 2. B6 2 /H 2 / 3 < u < 5p0t. 

By an adjustment of the constant C we can assume that u is a positive 
integer. Fix a density < p < 1 and an auxiliary density 

7/ 

X = P -. (5.27) 

H wet v ' 

Start with the basic coupling of three exclusion processes u > lu~ > rj with 
this initial set-up: 

(a) Initially {wj(0) : i ^ 0} are i.i.d. Bernoulli^) distributed and Wo(0) = 

1. 

(b) Initially w _ (0) = w(0) - 5 - 

(c) Initially variables {^,(0) : z 7^ 0} arc i.i.d. Bernoulli(A) and r/ (0) = 0. 
The coupling of the initial occupations is such that Wj(0) > 77^ (0) for all 
i ^ 0. 

Recall that basic coupling meant that these processes obey common Pois- 
son clocks. 

Let Q{t) be the position of the single second class particle between uj(t) 
and cj~(t), initially at the origin. Let {Xi(t) : i € Z} be the positions of 
the w — T] second class particles, initially labeled so that 

• • • < X_ 2 (0) < X_i(0) < X o (0) = < Xi(0) < X 2 (0) <■■■ 
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These second class particles preserve their labels in the dynamics and stay 
ordered. Thus the u(t) configuration consists of first class particles (the 
rj(t) process) and second class particles (the Aj(i)'s). P denotes the prob- 
ability measure under which all these coupled processes live. Note that the 
marginal distribution of (lu, u~ , Q) under P is exactly as it would be under 
P^. 

For x G Z, J"(i) is the net current in the w-process between space- 
time positions (1/2,0) and (x + 1/2, t). Similarly J£(t) in the ^-process, 
and Jx~ v (t) is the net current of second class particles. Current in the 
w-proccss is a sum of the first class particle current and the second class 
particle current: 

■£(*) = ■#(<) + ■£-"(*). (5.28) 

Q(t) <G {Xj(t)} for all time because the basic coupling preserves the 
ordering uJ~(t) > T)(t). Define the label mq(t) by Q(t) = X mQ ^(t) with 
initial value niQ(O) = 0. 

Lemma 5.4. For all t > and k > 0, 

P{m Q (t) >k}< e- 26k . 



Proof of Lemma 5.4- In the basic coupling the label rriQ(t) evolves as fol- 
lows. When X mQ -i is adjacent to X mq , niQ jumps down by one at rate 
p. And when X mq+ i is adjacent to X mQ , uiq jumps up by one at rate 
q. When X mQ has no X-particle in either neighboring site, the label tuq 
cannot jump. Thus the situation is like that in Lemma 5.2 (with a re- 
versal of lattice directions) with environment given by the adjacency of 
X-particlcs: u(m,t) — l{X m (t) = X m _\(t) + 1}. However, the basic cou- 
pling mixes together the evolution of the environment and the walk tuq, 
so the environment is not specified in advance as required by Lemma 5.2. 

We can get around this difficulty by imagining an alternative but dis- 
tributionally equivalent construction for the joint process (r),u)~, ui). Let 
(rj, u) obey basic coupling with the given Poisson clocks {N x ^ x±1 } at- 
tached to directed edges {x, x ± 1). Divide the u> — rj particles further into 
class II consisting of the particles us~ — rj and class III that consists only 
of the single particle oj — ui~ = 5q. Let class II have priority over class 
III. Introduce another independent set of Poisson clocks {N x ^ x±1 }, also 
attached to directed edges [x, x ± 1) of the space Z where particles move. 
Let clocks {N x ~* x±1 } govern the exchanges between classes II and III. In 
other words, for each edge {x,x + 1} clocks N x ^ x + 1 and N x + 1 ~ > - X are ob- 
served if sites {x,x + 1} are both occupied by uj — i] particles. All other 
jumps arc prompted by the original clocks. 

The rates for individual jumps are the same in this alternative construc- 
tion as in the earlier one where all processes were together in basic coupling. 
Thus the same distribution for the process (?y,cj~,cj) is created. 
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To apply Lemma 5.2 perform the construction in two steps. First con- 
struct the process (77, uj) for all time. This determines the environment 
u(m,t) = l{X m (t) = X m -i(t) + 1}. Then run the dynamics of classes II 
and III in this environment. Now Lemma 5.2 gives the bound for ttiq. □ 



Let u be a positive integer and 

,,2 



JZQOt 

By assuming t > C(B)8~ 1 we guarantee that 



> 1 and — — > > 4. 

A09t 40 



(5.29) 



Then 



k > — — > 4. (5.30) 



We begin a series of inequalities. 

P{Q(t) > V p t + u} 

< P{m Q (t) >k} + P{ J? VPti+u (t) - J L Vnj+„(*) > - fc >- ( 5 - 31 ) 

To explain the inequality above, if Q(t) > V p t + u and m,Q(t) < k then 
X k (t) > lV p t\ + u. This puts the bound 

^[VftJ+uW > 

on the second class particle current, because at most particles X\, . . . , X^-i 
could have made a negative contribution to this current. 

Lemma 5.4 takes care of the first probability on line (5.31). We work on 
the second probability on line (5.31). 

Here is a simple observation that will be used repeatedly. Process w can 
be coupled with a stationary density-p process u/ p ' so that the coupled 
pair (w,u/ p )) has at most 1 discrepancy. In this coupling 

I^W~^ (P) (*)I < 1- (5-32) 

This way we can use computations for stationary processes at the expense 
of small errors. 

Recall that V p = H'{p). Let c\ below be a constant that absorbs the 
errors from using means of stationary processes and from ignoring integer 
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parts. It satisfies |ci| < 3. 



EJ^tj+uW-EJfV.tj+uW ( 5 - 33 ) 
= tH(p) - (H'(p)t + u)p - tH(X) + {H'{p)t + u)X + a 
= -ltH"{p)(p-\Y-u{p-\) + Cl 
= t6(p - A) 2 - u(p - A) + d 
= t6(p - A) 2 - u{p - A) + ci + k - k 

2 2 2 

< k 

~ 100*6* ~ 10*0 + 20*0 ~ 

„.2 

fc. (5.34) 



25*6* 



The —3 in the definition (5.29) of k absorbed c\ above. 

Let X = X — EX denote a centered random variable. Continuing with 
the second probability from line (5.31): 



2 

< P|7jVftJ+«(*) _ JlVetj+uit) > 25^} 

< ^! V ar{j ryptJ+u (t) - J? VPtl+u (t)} 

< ^(Var{ J^ tJ+ Ji)} + Var{J^ J+ J*)}). (5.35) 

C is a numerical constant that can change from line to line but is indepen- 
dent of all the parameters. 

We develop bounds on the variances above, first for J u . Pass to the 
stationary density-p process via (5.32) and apply (5.11): 

Var{j^ PfJ+tt (*)} < 2V^{j [VPti+u (t)} +2 

= 2/9(1 - p)E\Q(t) - [VH\ -u\+2 

<E|Q(*) -V p t\ + u + 3 

<*(*) + 4u. (5.36) 

Let Var A denote variance in the stationary density-A process and let 
Q n (t) denote the position of a second class particle added to a process r/. 

Var{jf v , tJ+lt (t)} < 2Var A {J LyPtJ+n (i)} +2 

< E A |Q''(*) - [V p t\ -u\ + 2 

< E A |Q'"(*) - V p t\ +4u 
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Introduce process (£~(t), Q^(t), r](t), Q v (t))t>o coupled as in Theorem 5.6, 
where £ starts with Bernoulli(p) occupations away from the origin and 
initially Q^(0) = Q v (0) = 0. Below apply the triangle inequality and use 
inequality Q^(t) < Q v (t) from Theorem 5.6. Thus continuing from above: 

= E|Q"(i) - Q c (t) + Qt(t) -Vt\+ 4u 
< E{Q r '(t)-Qt(t)}+-E\Q((t)-V p t\+4u 
= V x t-V p t + ^(t)+4u 
= 26t{p - A) + + 4u 

= *(i) + 5«. (5.37) 

Marginally the process (C, Q^) is the same as the process (ui,Q) in the 
coupling of this section, hence the appearance of ^(t) above. Then we 
used (5.10) for the expectations of the second class particles and the choice 
(5.27) of A. 

Insert bounds (5.36) and (5.37) into (5.35) to get 

P{^ J+ .(«) ^ J[vn^)> -k}<Ce\^{t) + ^)- (5-38) 
Insert (5.30) and (5.38) into line (5.31) to get 

P{Q(t) < Vt ~u}< CO 2 ( L*( t ) + ^)+ e -« 2 /20t (5 39) 

and we have verified (5.24) for Case 2. 

Combining (5.39) and (5.26) gives the conclusion of Lemma 5.3. □ 

Next we extend the bound to both tails. 

Lemma 5.5. Let B € (0,oo). Then there exists a numerical constant 
C G (0,oo) and another constant cq{B) £ (0,oo) such that, for all densities 
< p < 1 andt> cq(B)6~ 1 , 

P"{ \Q(t) - V»t\ > u} 

\C9 2 ( Lm) + ^)+ 2e-" 2 /ct ; B9 2 /H 2 ^ <u< 20t/3 (5.40) 
< J V u u A J 

\2e- u / c , u>20t./3. 



Proof. The corresponding lower tail bound is obtained from (5.24) by a 
particle-hole interchange followed by a reflection of the lattice. For details 
we refer to Lemma 5.3 in [BS09a]. □ 
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Proof of the upper bound of Theorem 5.4- We integrate (5.40) to get the 
bound (5.13) on the moments of the second class particle. First for m = 1. 

/>oo 

*(t) = / P"{ \Q(t) - Vt\ > u} du 
Jo 

roc .2 j-2 . 

< B6 2 ' 3 t 2 ' 3 + C9 2 ( + ^r)du 

J Be 2/3 t 2/3 \M 4 U A J 



e~ u2/ct du + 2 / e~ u / c du 

9 2 /3 t 2/3 J20t/3 



DO 



< — -^(t) + {b + — -)e 2/3 t 2 

- 3B 3 w V 2B 2 / 

+ 02/3 e +ZL,e 

Ci(-B) is a new constant that depends on Set B = C 1//3 to turn the 
above inequality into 

*(*) < — « /3 * /3 + -^73- ex P( c x ) + 

The second term on the right above forces us to restrict t further. We can 
fix a constant Co large enough so that, for a new constant C, 

< C6 2/3 t 2/3 provided i > c 6»^ 4 . (5.41) 

Restrict to t that satisfy this requirement and substitute this bound on 
\I/(i) into (5.40). Then upon using u > B0 2 / 3 t 2 / 3 and redefining C once 
more, we have for B6 2 / 3 t 2 / 3 <u< 20i/3: 

P"{ \Q(t) - VH\ >u}< C^-Ar + 2e~ u2/ct . (5.42) 

Now take 1 < m < 3 and use (5.42) together with the second case of 
(5.40): 

E p |Q(t) - V p t\ m =m P p { \Q(t) - Vt\ > uju™- 1 du 
Jo 

/>oo 

< B m e 2m ' 3 t 2m ' 3 + CmQ 2 t 2 I u™- 4 du 

92/3^2/3 



OO 



+ 2m / e~ u ICt u m - x du + 2m e" u/c u" 1 " 1 du. 

Jse 2 /3t 2 /3 J20t/3 

Performing and approximating the integrals gives 

e p |QW - vnr < -^—e 2m ' 3 t 2m ' 3 

3 — m 

provided t > co# -4 for a large enough constant cq. □ 
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5.6 Proof of the lower bound for second class particle 
moments 

By Jensen's inequality it suffices to prove the lower bound for m = 1. 
Let Cub denote the constant in the upper bound statement that we just 
proved. We can also assume cq > 1. Fix a constant b > and set 

ai=2C UB + l and a 2 = 8 + V326 + 8 \JC UB ■ 

Increase b if necessary so that 

b 2 - 2a 2 > 1. (5.43) 

Fix a density p G (0, 1) and define an auxiliary density A = p— bt~ 1 ' 3 9~ 1 ' s . 
Define positive integers 

u = [ ai t 2/3 2/3 \ and n = [V x t\ - \ V p t\ + u. (5.44) 

By taking cq large enough in the statement of Theorem 5.4 we can ensure 
that A G (p/2, p) and u G N. 

Construct a basic coupling of three processes 77 < r/ + < £ with the 
following initial state: 

(a) Initially 77 has i.i.d. Bernoulli(A) occupations {?7;(0) : i =/= —n} and 

»?-n(0) = 0. 

(b) Initially r] + (0) = rj(0) + S- n . Q^ nS> {t) is the location of the unique 
discrepancy between rj(t) and ?7 + (i). 

(c) Initially ( has independent occupation variables, coupled with 77(0) 
as follows: 

(c.I) 0(0) = r)i(0) for -n < i < 0. 
(c.2) C-n(0) = 1. 

(c.3) For i < —n and i > variables Ci(0) are i.i.d- Bernoulli(/3) and 

G(o)>/?i(o). 

Thus the initial density of £ is piecewise constant: on the segment {— n + 
1,...,0} C(0) is i.i.d. with density A, at site — n £(0) has density 1, and 
elsewhere on Z f(0) is i.i.d. with density p. The reason for the gap in the 
C — 77 second class particles across (— n, 0] is to get an upper bound on the 
second-class particle current that is not too large for subsequent arguments 
((5.48) below). 

Label the ( — 1] second class particles as {Y m (t) : m G Z} so that initially 

• • ■ < y_i(o) < y (o) = -n - Q ( - n) (o) < < y x (0) < f 2 (o) < ■ ■ ■ 

Let again rag (rj) be the label such that Q( _n) (rj) = F„ lQ(t) (t). Initially 
mg(0) = 0. The inclusion Q^~ n \t) G {Y m {t)} persists for all time because 
the basic coupling preserves the ordering £(f) > f] + (t). Through the basic 
coupling rriQ jumps to the left with rate q and to the right with rate p, 
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(5.47) 



but only when there is a y-particle adjacent to Y mq . As in the proof of 
Lemma 5.4 we can apply Lemma 5.2 to prove this statement: 

P{m Q (() < -k} < e- 28k . (5.45) 

By the upper bound already proved and by the choice of a\ , 

P{Q(-")(i) > \ V p t\} = P{Q(-»)(t) > —n + [V x t\ +u] 

< u^\Q^\t) - n - [V*t\ | < (5-46) 
For the complementary event we get a lower bound: 

§ <P{Q(- n >(*) < \yn\] 
< P{m(t) < -k} + P{ jf VPt] (t) - J* VPt] (t) < k}. 

The reasoning behind the second inequality above is this. If Q^~ n '(t) < 
[Vt\ and m Q (t) > -k then Y-k(t) < [V p t\. This implies a bound on the 
second class particle current: 

J [vn\ (*) ~ J lvn\ W = J fv»'t} (*) ^ ( 5 - 48 ) 
Put k = La 2 t 1/3 1/3 J - 2. Then by t > 6~ A and the definition of a 2 , 

P{m Q {t) < -k} < e" 2 < 1/4. (5.49) 
Combine (5.47) and (5.49) and split the probability: 

\ < P{J C lVHl (*) - (*) ^ «2 <1/30l/3 - 2 } 

< P{ J c [VPt] (t) < 2a 2 t l ' 3 d 1 ' 3 + t 0( 2p A - A 2 )} 

+ P{ J^ pfJ (t) > a 2 t 1/3 # 1/3 + *0(2pA - A 2 ) + 2}. (5.50) 

Consider next line (5.50). The 77-process can be coupled with a station- 
ary P A -process with at most one discrepancy. The mean current in the 
stationary process is 

E x [J lVPtl (t)} = tH(X)-X[VH\ 

< tH(X) - XV p t + 1 = t6(2 P X - A 2 ) + 1. 

Hence 

line (5.50) < P x { Jy VPt] (t) > a 2 t 1/3 6 1 / 3 + t6{2pX - A 2 ) + 1} 
< P x {l [VPti (t) > oatVafli/a} < a^ 2 t- 2 ^e- 2 ^V a r x [j [VPti (t)} 
E x \Q{t) ~ [VPt\\ E x \Q(t) — VH\ , 2b t 1 



a 2 t 2 / 3 9 2 / 3 a 2 t 2 / 3 9 2 / 3 a\ a 2 t 2 l 3 B 2 l 3 

<C UB a^ + ^ + ^<\. (5.51) 
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After Chebyshev above we applied the basic identity (5.11) for which we 
introduced a second class particle Q(t) in a density-A system under the 
measure P A . Then we replaced [V p t\ with V x t and applied the upper 
bound and properties of a 2 . 

Put this last bound back into line (5.50) to be left with 

\ < P{J C [VPti (t) < 2a 2 t 1/3 d 1 / 3 + te(2p\ - A 2 )}. (5.52) 

Next we replace the (^-process with a stationary density- p process by 
inserting the Radon-Nikodym factor. Let 7 denote the distribution of the 
initial £(0) configuration described by (cl)-(c3) in the beginning of this 
section. As before v p is the density-p i.i.d. Bernoulli measure. The Radon- 
Nikodym derivative is 

/(w)=l£(w)=ii{w_ n =i}- n (-ik=i}+^iv.=o})- 

V " i=-n+l ^ " 

Bound its second moment: 

^(/ 2 ) = ifl + %P^y<p- 1 e^) 2 Mi-rt < C2{p) (5.53) 
p\ p{l-p)J 

where condition t > cq9~ 4 implies a bound C2(p) < 00 independent of t 
and 9. 

Let A denote the exclusion process event 

A = {Jyvn\ (*) < 2a 2 t 1 ' 3 d 1 / 3 + t0{2p\ - A 2 )}. 
Then from (5.52) 

§ < P{C G A] = J P"{A) 7 (dw) = J P u (A)f{Lu) v»(du>) 
< {P"{A)) 1/ \E"{f)f /2 < c 2 (p) 1 / 2 (P"(A)) 1/2 . (5.54) 
Note the stationary mean 

E p [ J [VPt] (<)] = tH(p) - p[V"t\ = tBp 2 + pV p t - p[V p t\ > tBp 2 . 
Continue from line (5.54), recalling (5.43): 

(64c 2 (p))- 1 < P P (A) = P p {Jy V pt\ (t) < 2a 2 t 1 ' 3 6 1 ' 3 + t£(2p\ - A 2 )} 

< P p { J lVPtl (i) < 2a 2 t 1 ' 3 B 1 ' 3 - t9{p - A) 2 } 
= P p { ~J\yn\ (t) <-(b 2 - 2a 2 )t 1 ' 3 9 1 l 3 } 

< t -a/3g-a/3 Var^[J LV ptj(*)] < r 2/3 <T 2/3 E"|Q(t) - Vt\. 
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This completes the proof of the lower bound. We finish with some ob- 
servations about the need for the two key assumptions, asymmetry and 
H"{p) ± 0. 

For symmetric SEP 9 = and consequently the Chebyshev step above 
cannot be taken. 

To observe where H"(p) < came in we need to backtrack a little. At 
stage (5.52) we have the inequality (ignoring now small errors due to integer 
parts etc.) 

| < P{-/fv-tj(*) ^ 2a 2 t 1 / s 6 1 / 3 + E x (J lVPti (t))}. 

The Radon-Nikodym and Schwarz trick turned this into an inequality for 
a stationary process: 

< c < P p {Jy V n\ (t) < 2a 2 t 1/3 e 1/3 + E x (J [VPt] (t))} 

= P p {~Jyvn\ (t) < 2a 2 t 1 ^9 1 ^ + E x (J lVPt} (i)) - EP(J lVHi (*))}. 

(5.55) 

Compute the means on the right-hand side inside the probability, remem- 
bering that V p = H'(p) and Taylor expanding H(X): 

E x (J [VPti (t)) - EP(J lVPti (t)) = t[H(X) - XH'ffi) - H{p) + P H'{p)\ 

= t[\H"{p){\-pf + o{e\\-p\ 3 )} 

= -03^t 1/3 fl 1/3 + 0(1) 

with ^H"(p) = —a 3 9 < in the last step. Thus the constants can be 
adjusted so that the probability in (5.55) is a deviation. Chebyshev can be 
applied to conclude that the current variance is of order i 2 / 3 6» 2 / 3 . But if 
H"{p) = there is no deviation to take advantage of. 

Further comments and references 

The proofs of this chapter are based on [BS09a] . This article gives simpler 
proofs for the results in [BS07] and [BS09b]. Precursors of these variance 
bounds were first proved for last-passage models that correspond to totally 
asymmetric versions of particle systems: in [CG06] for the Hammersley pro- 
cess and in [BCS06] for the corner growth model associated with TASEP. 

The Tracy- Widom type limit distribution for TASEP current was first 
proved for the step initial condition in [ JohOO] , then for the stationary case 
(Theorem 5.5) in [FS06]. The larger picture of TASEP fluctuations from 
various initial conditions is presented in [BAC09]. 

Another line of work has produced comparison theorems that allow one 
to conclude that Laplace transforms of t~ 1 'Var p [Q(t)} for different asym- 
metric exclusion processes are within constant multiples of each other. In 
this sense, for the order of this Laplace transform there is universality for 
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all finite range asymmetric exclusion processes. These results come from 
resolvent techniques [Set03, QV07, QV08]. 

The central limit theorem for the current in directions other than the 
characteristic V was proved first by Ferrari and Fontes [FF94]. This was 
generalized to other particle systems such as certain zero-range and brick- 
layer processes by Balazs [Bal03] . 

Consider symmetric simple exclusion, namely the case p = q = 1/2. 
Then V p = 0. Equation (5.11) together with the observation that the sec- 
ond class particle is a simple symmetric random walk tell us that Var p [Jo(i)] 
is of order i 1 / 2 , exactly as for independent particles in Section 2. And in- 
deed the current process does converge to fractional Brownian motion (see 
[PS08] and its references). 

6 Zero range process 
6.1 Model and results 

From the perspective of universality it would be highly desirable to extend 
the results of Section 5 beyond exclusion processes. Throughout the 40- 
year history of the subject of interacting particle systems, the zero range 
process has been a much-studied relative of the exclusion process. In this 
section we indicate how the bounds for second class particles and current 
variance are proved for a class of totally asymmetric zero range processes 
(TAZRP) with concave jump rate functions. 

Definition and graphical construction. In contrast with the exclu- 
sion process, the zero range process does not restrict the number of particles 
allowed at a site. The state of the process at time t is r)(t) = (?7i(i))iez £ ^+ 
where Tji(t) € Z + denotes the number of particles present at site i at time 
t. We consider the case where particles take only nearest-neighbor jumps 
to the right. 

Each zero range process is characterized by a jump rate function g : 
Z + — > M + . It automatically has the value g(0) = 0. The meaning of g for 
the process is that when the current state is n, g(j]i) is the rate at which 
one particle is moved from site i to i + 1. You can interpret this as saying 
that each of the particles at site i jumps independently with rate rj~ l g(rji) 
or that some particular one moves next (say, the bottom one is moved to 
the top of the next pile) at rate g(r]i). It is immaterial for we do not label 
our particles. (Except again we will label certain second class particles as 
we did for ASEP, but we will come to that later.) These jump events take 
place independently at all sites, exactly as for ASEP. 

We shall assume that 

g is nondecreasing and < g(k) < 1 for k > 0. (6.1) 
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With a bounded g we can perform the following concrete construction of 
the process rj(t) in terms of independent rate 1 Poisson processes {Ni} and 
i.i.d. Uniform(0, 1) variables {U^k}- Attach clock N to site i, and give each 
jump time of Ni its own U^. Now if t is a jump time for JV* with its uniform 
Ui k, then move one particle from i to i + 1 if Ui k < g(T]i(t— )), otherwise 
not. Repeat this step at all sites and all jump times. The result is the 
desired one: independently at each site i, a jump from i to i + 1 occurs in 
the next small time interval (t, t + dt) with probability g(r)i(t))dt + 0(dt 2 ). 
The generator of this TAZRP is 



that acts on bounded cylinder functions tp on Z+ and rf' l+1 = r\ — Si + 
Si + i. We will not use the generator in the text. It can be used to check 
the invariance of certain distributions on the state space, as for ASEP in 
Remark 5.2. 

Invariant distributions. Part of the reason for the popularity of ZRP 
is that, just like ASEP, it has i.i.d. invariant distributions. We denote these 
by {^ p }o< P <oo indexed by density p = E p {r]i). 

Here is the definition of these measures. Let 9 denote a real parameter, 
and on Z + define a probability distribution 

1 p ek 



Zg g{k)\ 



defined for 9 such that 



e 6k 



, *(*)! 



< CO. 



Here <?(0)! = 1 and g(k)\ = g(l) ■ ■ ■ g(k) for k > 0. Define the mean density 
function p(9) = kX e (k). It is smooth and strictly increasing on the 
open interval where Zg < oo. Let its inverse function be 9{p) and then 
reparametrize the distributions in terms of density: 

1 e e ^ k 
Ze( P ) 9{k)\ 

Finally, the actual invariant measures for ZRP are the product measures 
on the state space Z5.: 

V P (dTj)=6$ V S(dT H ). (6.2) 



We write P p for probabilities and E p for expectations for the stationary 
process whose marginal rj(t) has distribution v p . 
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Basic coupling and second class particles. Basic coupling works 
exactly as it did for exclusion processes: two (or more) zero range processes 
obey a common set of Poisson clocks {-/Vj} and uniform variables 

We write a boldface P for the probability measure when more than 
one process are coupled together. In particular, P p represents the situ- 
ation where the initial occupation variables 77^ (0) = 77^" (0) are i.i.d. mean-p 
Bernoulli for i 7^ 0, and the second class particle Q starts at Q(0) = 0. 

More generally if two processes 77 and to are in basic coupling and w(0) > 
?7(0) (by which we mean coordinatewise ordering Wi(0) > 77* (0) for all i) 
then the ordering ui(t) > 77(f) holds for all < f < 00. The effect of 
the basic coupling is to give priority to the 77 particles over the uj — 77 
particles. Consequently we can think of the w-process as consisting of 
first class particles (the 77 particles) and second class particles (the u — 77 
particles) . 

Current. The current is defined as for ASEP: for x € Z and f > 0, J x (f) 
is the net left-to-right particle current across the straight-line space-time 
path from (1/2,0) to (x+1/2,*). 

The flux function is again 

H(p) = _E p [rate of particle flow across a fixed edge] = E p [g(rji)]. 

Expectations of currents can be computed as for ASEP: 

E p [J x {t)] = tH(p) -xp, x e Z, f > 0. (6.3) 

The characteristic speed at density p is defined the same way as before: 

V p = H'{p). (6.4) 

As for ASEP, the first task is to establish the identities that connect 
current variance and the second class particle. These identities for ZRP 
develop the same way as for ASEP, except that a new initial distribution 
for the coupled process appears. Define a probability distribution on 
Z + by 

°o( k ) = VarP( ] Yl (m-pH{m), fceZ+. 

V IQ) m=k+l 

Define a product distribution v p on the state space Z^. that obeys the 
marginals v§ of the stationary distribution at all sites except at the origin 
where the distribution is : 

v p {dr]) = ((&v p {dr ]i ))®l> p {dr 1o ). 

Let P p be the probability distribution of a pair (77, 77+) that satisfies 77+ (f) = 
77(f) + Sqm (so there is one discrepancy), obeys basic coupling, and whose 
initial distribution is such that 77(0) ~ v p and 77+ (0) = 77(0) + So (in other 
words, Q(0) = 0). 
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Theorem 6.1. For any density 0<p<oo, z£Z and t > we have these 
formulas. 

V a r»[J z (t)} = J2\3 - A Cov^feCt), 770(0)], (6.5) 
Cov"[ Vj (t), rj o (0)] = v^( Vo )P p {Q(t) = j}, (6.6) 

and 

W[Q(t)] = V p t. (6.7) 

Equation (6.5) is proved the same way as for ASEP. Equation (6.6) and 
the definition of z/q come from a short calculation which we show below in 
Section 6.2. We omit the proof of (6.7). Formulas (6.5) and (6.6) combine 
to give the key equation that links current variance with the second class 
particle: 

Var p [J 2 (t)] = Var p (r/ )E>|Q(t) - z\. (6.8) 

Next we state the main result which again consists of upper and lower 
moment bounds for a second class particle, this time under the measure 
P p . We need a significant restriction on the concavity of the jump rate g: 

3 < r < 1 such that g(k + 1) - g(k) < r(g(k) - g{k - 1)). (6.9) 

A class of examples satisfying this hypotheses is given by g(k) = 1 — 
cxp(— ak b ) with a > 0, b > 1. g(k) can also be constant from some fco 
onwards. 

Theorem 6.2. Fix a density < p < oo and consider a pair of coupled 
ZRP's under the measure P p . Assume the jump rate function g satisfies 
(6.1) and (6.9). Then for 1 < m < 3, large enough t € R+ and a constant 
C, 

1 t 2m/3 < gp r ,Q/ t s _ ypnm j < _^_ t 2m/3_ /g 1Q ) 

(7 3 — m 

The constants C and how large i needs to be may depend on the density 
p. Combining (6.5) and (6.10) gives the bounds for the variance of the 
current seen by an observer traveling at the characteristic speed V p : for 
large enough t, 

c -i t 2/3 < V ar p [J L y Pt j(i)] < Ct 2/3 . 

For the remainder of this chapter we discuss parts of the proof. Once we 
have the fundamental identities that tic together moments of the second 
class particle and the variance of the current, the proofs for the upper and 
lower bounds given for ASEP in Sections 5.5 and 5.6 can be adjusted to 
work for TAZRP. We will not repeat those derivations. Instead, we focus 
on the key ingredients that made the argument work for ASEP, and discuss 
how to provide these ingredients for TAZRP. There are two key points that 
we need in order to repeat the argument for TAZRP: 
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1. We need a construction that includes a given second class particle as 
a labeled member of a density of second class particles, and then we 
need a tail bound for the label as the one given for ASEP in Lemma 
5.4 and (5.45). 



2. We need a coupling that keeps the second class particle of a system 
with higher density behind the second class particle of a system with 
lower density. For ASEP this was Theorem 5.6, which was used to 
obtain (5.37) for the ASEP upper bound. 



We turn to these points in Section 6.3 below, after developing the variance 
formula. 



6.2 Variance identity 

Define F(-l) = and 

Interpret below Vq(— 1) as 0. 

Cov P Mt), m (0)} = E"Mt)(r, (0) -P)] = J2 E h(t)\vo(0) = k](k- p)v?(k) 

k>0 

= £ E[ Vi (t)\r, (0) = k](F(k - lK(fc - 1) - F(kX(k)) 

k>0 

= J2Eh(t)\vo(0) =k + l]F{kX(k) -Y, E h(t)\vo(0) = k}F(kX(k). 

k>0 k>0 



Construct a coupling of rj + (t) = r](t) + Squ) with the discrepancy initially 
at the origin Q(0) = 0, and so that r/(0) has ^-distribution. Then, due to 
the product form of the initial distribution, 



E[77, (t)\vo(0) = k + l}= E[77+(i)|77+(0) = k + 1] = E[rj+ (t)\ m (0) = k]. 
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Then continuing from above, 

Cov"^), 770(0)] = ^EK+^MO) = k]F(kX(k) 

k>0 

-Y,nm(t)\M0) = k]F(kX(k) 

k>0 

= 5>fo+(t) m(t)\vo(o) = k]F(kK(k) 

k>0 

= E[(r ) +(t)-r H (t))F(r 10 (0))] 

= E[l{Q(t) = i}F(r, o (0))] 

= J2nHQ(t) = i}\m(0) = k]F(k)v p Q {k) 

k>0 

= Var'M £ E[l{Q(t) = i}|r/ (0) = k]v p {k) 

k>0 

= VaT p ( Vo )P p {Q(t) = »}. 

This proves (6.6). 

6.3 Coupling for the zero range process 

Next we describe a coupling of two processes with labeled discrepancies 
(second class particles) between them, and then two randomly evolving 
labels that achieve simultaneously both goals (1) and (2) mentioned above. 
This construction works for any TAZRP with concave jump rate function 
g. Getting tail bounds for the label processes is the serious bottleneck of 
this proof, and that is where we need the restrictive assumption (6.9). 

Let two processes rj < uj evolve in basic coupling. This pair (77, ui) to- 
gether with the labeled and ordered lo — r\ second class particles ■ • • < 
X- 2 (t) < X_!(t) < X (t) < Xx(t) < X 2 (t) < ■ ■■ form a "background" 
process on which we define two label processes y(t) and z(t). The lu — i] 
second class particles are kept in order by requiring that, whenever a sec- 
ond class particle jumps to the right, the X-particle with highest label is 
moved. 

The label processes will satisfy y(t) < z(t), we will be able to bound y{t) 
stochastically from above, z{t) stochastically from below, and the following 
two pairs of processes will individually be in basic coupling: 

(u>~, to) = (u> - S Xy , u) and (??,?7 + ) = (77,77 + S X J- (6-H) 

The definition of the label processes is partly forced on us by the require- 
ment that jumps of X y and X z must replicate the rates required by the 
basic coupling. Additionally, we devise the joint process (y, z) so that the 
order y < z is maintained. 
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The rule is that after each jump among {oj,rf) that in any way affects 
the site where X y resides, the value of y is refreshed randomly. Let a and 
b denote the minimal and maximal labels at the site i where X y resides 
after the jump. If X y resides at a site other than X z , then y chooses a new 
value y' according to these probabilities: 



a with probability 



b with probability 



gjfJi - 1) -gfa) 
g(u>i) - g{m) 

g(ui)-g(rn) 



(6.12) 



If g(u>i) - g(rji) = then y' = a. 

Similarly, after a jump in the background process that affects the site 
where X z resides, if X z and X y are not together, then z takes the new 
value z' as follows (with b again the maximal label at the site i = X z after 
the jump): 



1 with probability 



with probability 



g(u)j) - fffa + 1) 
g(ui) - g(r)i) 

g(vi + 1) - gfa) 

g(ui) -g(r)i) 



b - ly 5(^1) - g(vd 

with probability 



If g(uJi) — g{r/i) = then z' = b. 

Finally, if after the jump X y = X z = i, then the labels are refreshed 
jointly as follows: 

a ^ with probability 9^)- 9 (Vi + l) 



bj g{ui) - g(r]i) 

If gfa) - g(r]i) = then (y',z') = (a, b). 

The jump rules preserve y < z. Note that marginally y' obeys proba- 
bilities (6.12), and similarly for z' . Concavity of g was used to define the 
middle case in the joint rule. 

Let us observe why these rules give the pair (lo~,u>) = (u> — Sx v ,oj) the 
same rates this pair would have in basic coupling. The requirement is 
that a jump across edge (i, i + 1) occur for both processes at rate g(ui~), 
and only for lo at rate g(u>i) — g(uj~). This requires thinking through a 
few cases. Only the site where X y resides needs attention since elsewhere 
, uj) jump together according to ZRP rates. 



bj g(u)i) - g(rji) 

_ g{uj) - g{ui - 1) 
g(wi) - g(r}i) 

^ with probability 9^) - g^ - 1) 
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(i) In the basic coupling of (t],uj), an (£, i + 1) jump occurs in r\ at rate 
g[j]i). Then both uj and uj~ experience this jump. 

(ii) An uj — rj second class particle jumps at rate g{uJi) ~ g{">]i)- Prior to 
this jump y chose the top label with probability given by the second 
line of (6.12), hence the rate at which X y jumps is 

(g(u)i) - g(r)i)) r— = g(uj t ) - g[ui - 1). 

9m) - gym) 

Thus at this rate uj experiences the jump but uj~ does not. 

If prior to this jump y chose the bottom label, then both uj and ui~ 
experience this jump, and this happens with rate 

i i s / \\ 9(ui - 1) -g(vi) , ^ / \ 

g[ui) - giVi) 

Adding up the rates we see that the rates of basic coupling have been 
realized. A similar argument is given for (77,77+) = (77,77 + Sx z )- 

We come to the unique point in the proof where assumption (6.9) is used, 
namely the tail bounds for the labels. 

Lemma 6.1. Let the labels start with y(0) = z(0) = 0. Under assumption 
(6.9) we have these bounds: P{y(t) > k} < r k and P{z(t) < — k} < r k for 
all k G Z + and t > 0. 

Proof. We do the proof for y(t). The bounds are valid conditionally on 
the evolution (77, uj) of the background process. So assume this background 
evolution given. Then we think of y(t) as an integer- valued Markov chain 
that is subject to jumps triggered by the background environment. Each 
jump happens at some site i with range of labels {a, . . . , b} and occupation 
variables uj% > 77, > that together satisfy 

— rji = b — a + 1. 

Given the current value y, the new value y' is obtained by the following 
rules, which of course are consistent with (6.12): if g(uji) — g(i]i) = then 

, \y if y<a or y>b 
ya if y e {a, . . . ,bj 

while if g(uji) —3(771) > then 

y if y < a or y > b 

a with probability — -—— if y G {a, . . . , b} 

g(uJi) - g{r)i) (6.14) 

b with probability if y e {a, . . . , b}. 

g(uji) ~ g{r}i) 
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There are infinitely many such jumps in any finite time interval, but all 
but finitely many leave y unchanged. To get around this difficulty, we can 
first freeze all the Poisson clocks outside space interval [— M, M], prove the 
lemma there, and then let M oo. Since rates are bounded, in any given 
bounded block of space-time the finite-M process agrees with the infinite 
process for all large enough M. 

To prove the lemma we show that every jump of type (6.13)-(6.14) pre- 
serves the geometric tail bound, regardless of the values a,b,uji,rji. So 
suppose y is an integer-valued random variable such that 

P(y >k)<r k for k > 0, 

and define y' via (6.13)-(6.14). We wish to show that P(y' > k) < r k for 
k > 0. 

The case (6.13) is clear since there y' < y. Let us consider the case 
g(uji) — g(r)i) > 0. Since the jump only redistributes the probability mass 
in {a, ...,&} to {a, b}, it suffices to check that 

P(y' >b)<r b (6.15) 

in the case b > 0. Using the jump rule (6.14), 

P(y , >b)=P(y' = b) + P(y'>b+l) 
g(ui) - g{u l - 1) 



g(u)i) - g(r)i) 



P(a < y < b) + P(y > b+ 1) 



If g(u>i) — g{uji — 1) = the conclusion (6.15) follows from the assumption 
on y. So we assume g(uii) — g(oJi — 1) > 0. Then by concavity all the 
^-increments between rji, . . . ,u>i are positive. Next write 

g(u)i) - g{r)i) ^ 

g{ui)-g{r]i) 

+ P{y > 6+1). 



For a < k < b 



e=uj,-b+k y v ; 3 w 

b g(uj, - b + k) - g(wi - b + k - 1) 

= {l-r)r — . 

g{uji) - g{uj t - 1) 
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Adding these up over a < k < b gives 
first term on the right in (6.16) 

gi^i) - g(vi) 

x ^ g(oJi — b + k) — g(cui — b + k — 1) 

^ g(wi) - gfai - i) 
Nfc ff(^i) - gfai - b + a - 1) „ 

= (1 — r)T ■ r ; r = (1 — TJT = T — T 

g(ui) - g(rji) 

Substitute this bound back up to (6.16) and use P(y > k) < r k twice: 
P(y' > b) < r b + P(y > b + 1) - r b+1 

g(Wi) - g{r)i) 

b g(u>j) - g(u>j - 1) ( a \ 

- r + — r~\ — 7 — s — [P{y >a)-r) 
g{ui)-g{r)i) 

<r b . 

Thus (6.15) has been checked and thereby the lemma has been proved for 
y(t). □ 
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